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{N| ' Abstract 

^V The Muskat problem involves nitration of two incompressible fluids throughout a 

porous medium. In this paper we shall discuss in 3-D the relevance of the Raylcigh- 
Taylor condition, and the topology of the initial interface, in order to prove its local 
existence in Sobolev spaces. 



1 Introduction 

0* . 

The Muskat problem (see ref. [13] and [2]) involves filtration of two incompressible fluids 
,£3 ■ throughout a porous medium characterized by a positive constant n quantifying its porosity 

and permeability. The two fluids, having respectively velocity fields v 3 , j = 1,2, occupy 
disjoint regions D 3 (D 2 = M 3 — D l ) with a common boundary (interface) given by the surface 
S = dD 3 . Naturally those domains change with time, D 3 = D 3 (t), as it does its interface 
S = S(t). We shall denote by p 3 (j = 1,2) the corresponding pressure and we will assume 
also that the dynamical viscosities p? and the densities p 3 are constants such that p} ^ /j, 2 , 

P 1 + P 2 - 
if} . The conservation of mass law in this setting is given by the equation V • v = (in the 

CO | distribution sense) where v = v 1 Xd 1 + v 2 Xd 2 - 

The momentum equation was obtained experimentally by Darcy |10|. [2] and reads as 

follows 

(*L v i = -Vpi- (0,0, pig), J = 1,2, 
K 

where g is the acceleration due to gravity. 

One can find in the literature several attempts to derive Darcy law from Navier-Stokes (see 
[TB] and [15] ) throughout the process of homogenization under the hypothesis of a periodic, 
or almost periodic, porosity. In any case the presence of the porous medium justify the 
elimination of the inertial terms in the motion, leaving friction (viscosity) and gravity as the 
only relevant forces, to which one has to add pressure as it appears in the formulation of 
Darcy's law. There are three scales involved in the analysis, namely: the macroscopic or bulk 
mass, the microscopic size of the fluid particle and the mesoscopic scale corresponding to the 
pores. In the references above one find descriptions of the velocity v as an average over the 
mesoscopic cells of the fluid particle velocities. Taking into account that each cell contains 
a solid part where the particle velocity vanishes, it is then natural to get the viscous forces 
associated to that average velocity, which is a scaled approximation of the laplacian term 
appearing in the Navier-Stokes equation. 

In this paper we shall consider the case of an homogeneous and isotropic porous material. 
Porosity is the fraction of the volume occupied by pore or void space. But it is important to 
distinguish between two kind of pore, one that form a continuous interconnected phase within 
the medium and the other consisting on isolated pores, because non interconnected pores can 



not contribute to fluid transport. Permeability is the term used to describe the conductivity 
of the porous media with respect to a newtonian fluid and it will depend upon the properties 
of the medium and the fluid. Darcy's law indicates such dependence allowing us to define 
the notion of specific permeability k and its appropriated units. In the case of anisotropic 
material k will be a symmetric and positive definite tensor, and then the methods of our proof 
can be modified to get local existence, but for a non homogeneous medium the properties of 
the tensor k(x) will have to be conveniently specified in order to have an interesting theory. 
The Muskat and related problems [H] have been recently studied O [161 El 13 E] • The 
first natural question asks for the evolution (existence) of such system, at least for a short 
time t > 0, and the persistence of smoothness of the interface S(t) if we begin with a smooth 
enough surface at time t = 0. One can deduce easily from this formulation that in the 
occurrence of such smooth evolution both pressures, modulo a constant, must coincide at the 
interface: 

P l \s(t) =P 2 \s(t)- 

Therefore we look at the case without surface tension (see article |11] where the regularizing 
effect of surface tension is considered). The normal component of the velocity fields must 
also agree at the free boundary: 

(v 1 - v 2 ) ■ u j = at S(t) 

where v 3 is the inner unit normal to S at the domain D 3 [y = —v 1 )- Furthermore the 
vorticity will be concentrated at the interface, having form 

curl(u) = ui(z)dS(z) 

where uj is tangent at S at the point z and dS{z) is surface measure. 

The main purpose of this paper is to extend to the 3-dimensional case the results obtained 
in [5] for the case of 2 dimensions, namely proving local-existence in the scale of Sobolev spaces 
of the initial value problem if the Rayleigh- Taylor condition (R-T) is initially satisfied (see 
[14] ) where this issue is studied from a physical point of view). In our case that condition 
amounts to the positivity of the function 

a = (Vp 2 - Vp 1 ) ■ (y 2 - v l ) 

at the interface S. Let us indicate that the R-T property also appears in other fluid interface 
problems such as water waves [6]. 

Together with that hypothesis, one also assume that the initial surface S is connected and 
simply connected. In the presence of a global parametrization X : M. 2 — > S, the preservation 
of that character will be controlled by the gauge 

F ^'^> = \X(a)-X W l|FWI1 - = gg |*M-*(ffl| K °°- 

Section 2 of this paper contains the deduction of the evolution equations for the interface 
S. In section 3 we prove the existence of global isothermal parametrization as a consequence of 
the Koebe-Poincare uniformization theorem of Riemann surfaces in the geometric scenarios 
considered in our work, namely: either double periodicity in the horizontal variables or 
asymptotic flatness. Let us add that given the non-local character of the operator involved, 
to obtain a global isothermal parametrization is an important step in the proof, whose main 
components are sketched in section 4. 



In closing our system (section 2) we need to control the norm of the inverse operator 
(I+AD) -1 where T> is the double-layer potential and |A| < 1. It is well known from Fredholm 's 
theory that those operators are bounded on L 2 {S). However since the surface S = S(t) is 
moving, a precise control of its norm is needed in order to proceed with our proof. That is 
the purpose of section 5 where the estimates for the double-layer potential are revisited. 

In sections 6 and 7 we develop the energy estimates needed to conclude local-existence. 
Let us mention that at a crucial point (more precisely just at that step where the positivity 
of a(a,t) (R-T) plays its role), we use the pointwise estimate 9(x)A9(x) > ^A6 2 (x) of [I], 
with A = y/-A. 

In the strategy of our proof it is crucial to analyze the evolution of both quantities o and 
F (section 8) at the same time than the interface X and vorticity uj. There are several pub- 
lications (see pQ for example) where different authors have treated these problems assuming 
that the Rayleigh- Taylor condition is preserved during the evolution. Under such hypothesis 
the proof can be considerably simplified, specially if one also assume the appropriated bounds 
for the resolvent of the double layer potential respect to a moving domain, or the existence 
of global isothermal coordinates, etc... It is our purpose of going carefully over such items 
what is responsible for the more delicate and intricate parts of this paper. 

2 The contour equation 

We consider the following evolution problem for the active scalars p = p(x, t) and p = p(x, t), 
x G R 3 , and t > 0: 

Pt + v-Vp=0, 

Pt + v- Vp = 0, 
with a velocity v = (1)1,1)2,1)3) satisfying the momentum equation 

pv = -Vp- (0,0, p), (2.1) 

and the incompressibility condition V • v = 0, where, without loss of generality, we have 
prescribed the values k = g = 1. 
The vector (p, p) is defined by 

, S, ^_ J (/Ap 1 ), xeD\t) 

ip,p){X l ,X 2 ,X 3 ,t)-^ ( ^ 2)) , GZ) 2 (t)=i 3^1 (i)) 

where p 1 7^ p? , and p 1 7^ p 2 . Darcy's law (|2.ip implies that the fluid is irrotational in the 
interior of each domain D 3 and because of the jump of densities and viscosities on the free 
boundary, we may assume a velocity field such that 

curlv = oj(a,t)5(x — X(a,t)), 

where dD^(t) = {X(a,t) G R 3 : a = (cti,a 2 ) G R 2 }, i.e. 

<curli;,^>= / w(u,t) ■ (p(X(u,t))da, (2.2) 

for any <p : R 3 — > M 3 vector field in C^° 

The incompressibility hypothesis ( < V ■ V, ip >= — < v, V</? >= for any <p G C^°(1R 3 )) 
yields 

v 1 (X{a,t),t) ■ N(a,t) = v 2 (X{a,t),t) ■ N(a,t), 



with N(a,t) = d ai X(a,t) A d a2 X(a,t), and equation (|2.2|) gives us the identity 

uj(a,t) = (v 2 (X(a,t),t) -v 1 (X(a,t),t)) f\N(a,t). 

Defining the potential (p by v(x,t) = \7(ft(x,t) for x £ M 2 \ dD^(t), we get 

n(a,t) = 4> 2 (X(a,t),t) - 4> l (X(a,t),t), 

d ai Q(a,t) = (v 2 (X(a,t),t) - v 1 (X{a,t),t)) • d ai X, 
d a2 n(a,t) = (v 2 (X(a,t),t) - v 1 (X{a, t), t)) -d a2 X. 
Then, one has the equality 

uj(a, t) = (v 2 (X(a, t),t) - v 1 (X(a, t),t)) A (d ai X(a, t) A d a2 X(a, £)) 

and therefore 

uj(a,t) = d a2 Q(a,t)d ai X(a,t) — d ai Q(a, t)d a2 X(a,t), (2.3) 

implying that V • curl v = in a weak sense. 

Using the law of Biot-Savart we have for x not lying in the free surface (x ^ X(a,t)) the 
following expression for the velocity: 

„(*,*) = _jL f -^M-Aw^. 
V ; A7rJ R 2 \x-X(f3,t)\ 3 yH> H 

It follows that 

X t (a) = BR(X,u){a,t) + d(a)d ai X(a) + C 2 (a)d a2 X(a), (2.4) 

where BR is the well-known BirkhofF-Rott integral 

BR(X,„ ){a , t) = -±PV I ] mzm ?Aumi . (2.5) 

Next we will close the system using Darcy's law: 
Since 

Vcj) = v(x,t) - n(a,t)N(a,t)S(x - X(a,t)) 

we have 

< A(f>, tp >= - < V0, Vip >= / n(a, t)N(a, t) ■ V<p{X(a, t))da, 

taking <p(y) = — 1/(4tt|x — y\) one obtain <j) in terms of the double layer potential: 

m = ~ I { x ~^% ■ ^ («)«(«)**■ 

47ry R2 \x-X(a)\ 6 

Darcy's law yields 

Ap(x,t) = -div(p,(x,t)v(x,t)) -d X3 p(x,t), 

that is 

Ap(x,t) = P(a,t)5(x - X(a,t)), 

where P{a,t) is given by 

P(a, t) = (p 2 - p})v{X{a, t),t) ■ N(a, t) + (p 2 - p^N-^a, t), 



implying the continuity of the pressure on the free boundary. 

Next if x 7^ X(a,t), i.e. x is not placed at the interface, we can write Darcy's law in the 
form 

p<j)(x,t) = -p(x,t) - px 3 

and taking limits in both domains D J we get at S the equality 

(p 2 ^ 2 (X(a, t),t) - /iV (A>, t),t)) = -(p 2 - p 1 )X 3 (a, t). 
Then the formula for the double layer potential gives 



2 Att J&\X(a)-X(l3)\* 

that is 



where 



Q(a, t) - A^V(n)(a, t) = -2A p X 3 (a, t), (2.6) 

V(Q)(a) = ±PV f JW—gW • N(P)il(P)d0, (2.7) 

2vr J R2 |X(q)-X(/3)| j 

_ p 2 -p l _ p 2 -p 1 

A ^-^ + fx l and A P- fI 2 +jl l- 

And the evolution equation are then given by (|2.3j )-l l2.7|l . where the functions C\ and C 2 will 
be chosen in the next section. 

Furthermore, taking limits we get from Darcy's law the following two formulas: 

d ai n(a,t) + 2A li BR(X,u)(a,t) ■ d ai X(a,t) = -2A p d ai X 3 (a,t), (2.8) 

d a2 n(a,t) + 2A fM BR(X,uj)(a,t) ■ d a2 X(a,t) = -2A p d a2 X 3 (a, t). (2.9) 

3 Isothermal parameterization: choosing the tangential terms 

Although the normal component of the velocity vector field is the relevant one in the evolution 
of the interface, it is however very important to choose an adequate parameterization in 
order to uncover and handle properly the cancelations contained in the equations of motion. 
Fortunately for our task we can rely upon the ideas of H. Lewy [12], and many other authors, 
who discovered the convenience of using isothermal coordinates in different P.D.E. namely for 
understanding how a minimal surface leaves an obstacle, but also in several fluid mechanical 
problems. 

Let us recall that an isothermal parameterization must satisfy: 

\X ai (a,t)\ 2 = \X a2 (a,t)\ 2 , X ai (a,t) ■ X a2 (a,t) = 0, 



for t > 0. 

Next we define 



n ( s 1 f a l ~ ffl BR fa ' X fo ~ BR l3i ■ Xjh j g 
WW = 7T / I «I2 \Y 12 d > j 

2vr J R 2 \a- (3\ 2 \Xf3 2 \ 2 

1 f a 2 - /3 2 BRfo ■ Xfc + BR/3 2 ■ X^ 



(3.1) 



2tt7 m2 |«-/3|2 \XpJ2 



and 



If a2 — (h BR/3 2 ■ Xp 2 — BR^ ■ Xp x 



C2(a) = — r : 2 ^ t v , 2 n - dp 



02 



27T 7 M2 |«-/3P |X 

1 f a\ — fii BR^ ■ X/3 2 + BR/3 2 ■ Xp t 



(3.2) 
2vr7 R2 |a-/3|2 |X ft | 2 P ' 

That is X t = BR + C x X ax + C 2 X a2 and 

X ai t = BR ai + GiA aiai + C2^- Ql « 2 "I" Wai^-cti + ^2ai^a 2 i 
^■02* = BR a2 + CiA aiQ , 2 + C2^-a 2 02 + ^l«2^-oi "I" ^2Q2^a2- 

Denoting / = (IXqJ 2 — \X a2 \ 2 )/2 and g = X ai • X a2 we have 

/t = {BR ai -X ai — BR a2 -X a2 )+Cif ai +C2f a2 + (C2 ai — Ci a2 )g+2Ci ai f+(Ci ai — C2 a2 )|X Q2 | . 

The expressions for Ci and C2 yield the vanishing of the sum of the first and the last 
terms in the identity above. Therefore we get 

ft = C\f ai + C2Ja 2 + (C2ai — C\ a2 )g + 2C\ ai f. (3.3) 

Similarly we have 

9t = (BR a2 -X ai +BR ai -X a2 ) + Cig ai +C2ga2 + {Cla 1 +C2a 2 )9 — '2C2a 1 f + {Cl a2 +C2a 1 )\X ai \ , 

and 

9t = Cig ai + C 2 g a2 + {Ci ai + C 2Q2 )g - 2C 2ai f. (3.4) 



The linear character of equations (13.3|) and (13. 4p allows us to conclude that if there is a 
solution of the system Xt = BR + C\X ai + C2X a2 and we start with isothermal coordinates 
at time t = 0, then they will continue to be isothermal so long as the evolution equations 
provide us with a smooth enough interface. 

The fact that one can always prescribe such coordinates at time t = follows from the 
following argument: In the double periodic setting we have a C 2 simply connected surface, 
homeomorphic to the euclidean plane ]R 2 , which, by the Riemann-Koebe-Poincare uniformiza- 
tion theorem, is conformally equivalent to either the Riemann sphere, the plane or the unit 
disc. The sphere is easily eliminated by compacity, but we can also rule out the unit disc 
because the double periodicity assumption in the horizontal variables imply the existence of 
an abelian discrete subgroup of rank two in the group of conformal transformations, and that 
event can not happen in the case of the unit disc. 

Therefore we have an orientation preserving conformal (isothermal) equivalence 

(/) : R 2 — ► S. 

Since S is invariant under translations t v (x) = x + 2irv, v S 1? x {0} it follows that f v {z) = 
(f)~ l o t u o 4>(z) must be a diffeoholomorphism of C = M. 2 and, therefore, it has to be of the 
form 

fu(z) = a u z + b u , for certain a u ,b u G C. 

Clearly the family f v is generated by f x = /(i )0 ,o), /2 = /(o,i,o)- Let 

fi(z) = aiz + bi, f2(z) = a2Z + b 2 



we claim that a\ = 02 = 1. Suppose that |oi| < 1 then we get /" (z) = a\z + b\{\ + a,\ + .. + 
a™~ ) a sequence converging to j-^- contradicting the discrete character of the group action. 
On the other hand, if |ai| > 1 then since 

fi l i z ) = f{~i,o,o)(z) = - 

we get a contradiction with the sequence f^ n (z). Therefore we must have a\ = e 2m8 for some 
< 9 < 1. Assume that < 9 < 1 then 



r( n )(~\ - p 2mnS „ , u (-, , _2ttj0 , , 27ri(n-l)6K _ 27rm0 



/W(^) = e ™^ + 6i(l + e^ y + + e ^ n -W) = e 2mnti z + h=- 



1-e 



2irind 



e 2iri6 



That is the sequence f n (z) is bounded, \f n {z)\ < \z\ + J ^ e , and therefore it contains a 
converging subsequence contradicting again the discrete character of the action. That is, we 
must have f\{z) = z + b\ and, similarly, fo{z) = z + 62, allowing us to conclude easily the 
double periodicity of the isothermal parameterization <f>. 

In the asymptotically flat case we start with an orientable simply connected surface S so 
that outside a ball B in M 3 it becomes the graph of a C 2 -function X3 = y>{x\, X2) satisfying that 
\D a ip(x)\ = o(\x\~ ), for every N and \a\ < 2, in particular, the normal vector ^ = '== = 

u(x) is pointing out vertically , = = 3> ^ for |x| big enough. 

It is then well known that one can find isothermal coordinates whose first fundamental 
form X(a, f3)(da 2 + df3 2 ) converge asymptotically to the identity. 

Again by the uniformation theorem S must be conformally equivalent to either C or the 
unit disc. But since outside B the surface S is conformally equivalent to C — B {^\{x^ = 0} 
it cannot be also conformally equivalent to D — K, for any regular compact set K contained 
in the unit disc D, because the harmonic measure of the ideal boundary is 1 in the case of D 
and for R 2 . 

4 Outline of the proof. 

The proof of local existence requires the following: 

1) A connected and simply connected surface S = S(t) parameterized by isothermal 
coordinates 

X :R 2 — >R 3 , X = X(a,t) 

with normal vector N(a,t) = X ai A X a2 and gauge 

F(X)(a,0) = \P\/\X(a)-X(a-P)\, 

such ||F(X)||ioo < 00 and |||-ZV| -1 ||_l°o < 00. 

2) The positivity of 

a(a,t) = -{Vp 2 {X(a,t),t) - Vp 1 (A(a, t), t)) ■ N(a, t) 

= (fj 2 - fi^BRiX, u)(a, t) • N{a, t) + {p 2 - p x )N^{a, t), 

where the last equality is a consequence of Darcy's law after taking limits in both domains 
D 3 . This is the Rayleigh- Taylor condition to be imposed at time t = 0, being a part of the 
problem to prove that it remains true as time pass. 



3) The estimates on the norm of (I — XD) , |A| < 1, T> = double layer potential (section 
5) allows us to obtain the inequalities: 

\\n\\ Hk+1 < P(\\x\\ 2 k+1 + \\F(x)f LOO + IHivr 1 !!^), 
\\u\\ H » < P(\\x\\l +1 + \\F(x)f Leo + IHivr 1 !!^), 

for k > 3, where P is a polynomial function and the norm || • ||/% is given by 

\\X\\ k = \\X X - ai|| L3 + \\X 2 - a 2 \\ L 3 + \\X 3 \\ L 2 + ||V(X - (a, O))]^ , 

as in (|7.ip below, and || • \\ H j denotes the norm in Sobolev's space HK 

4) A control of the Birkhoff-Rott integral BR(X,u): 

\\BR(X,co)\\ Hk < CP(\\Xf k+1 + \\F(X)\\ho + IHiVr 1 !!^). 

for k > 3. 

5) Energy estimates: The properties of isothermal parameterizations help us to reorganize 
the terms in such a way that 

j t \\x\\ut) < p(\\x\\i(t) + \\F(x)\\ 2 Loa (t) + wim^w^it)) 



Et^T / r^^ d k ai X(a,t) ■ A(d k ai X)(a,t)da, 



where k > 4, |VA(a)| 3 = (\d ai X(a)\ 2 + \d a2 X(a)\ 2 ) 3 / 2 and A = (-A) 1 ^ = (R x (d ai ) + 
R2(d a2 ))- Then the pointwise inequality 

6A(6) - h(6 2 ) > 0, 

together with the condition a > allows us to get rid of the dangerous terms in the inequality 
above (i.e. those involving (k+l)-derivatives of X) to obtain the estimate 

j t \\X\\ 2 k (t) < P(\\X\\ 2 (t) + \\F(X)\\ 2 Lao (t) + liw^iu-ct)). 

6) Finally we need to control the evolution of ||F(A)||ioo(i) and inf(t) = inf a(a,t) 
which is obtained via the following estimates 

^||F(X)||i„(t) < P(\\X\\ 2 t (t) + ||F(X)||1„(() + IHWI-'lli-W) 

s ( J®' - « p(lwll(i) + » f w»-< 4 > + iim-iu-m). 

7) All those facts together yield the inequality 

f t E(t) < CP(E(t)), 
for the energy: 

E(t) = \\X\\ 2 k {t) + ||F(*)||£.(t) + IHiVr 1 !!^^) + (inf(t))- 1 
where k > 4, C is an universal constant and P has polynomial growth (depending upon k). 



At this point it is not difficult to prove existence of a solution, locally in time, so long as 
the initial data A(0) is in the appropriate Sobolev space of order k, k > 4, and the Rayleigh- 
Taylor and not-selfintersecting conditions (cj > c > 0, ||-F(X(0))||£°o < oo) are satisfied. 
Finally let us point out that since our existence proof is based upon energy inequalities an 
extra argument is needed to prove uniqueness. Nevertheless that task is much easier than 
proving existence ( the interested reader may consult the forthcoming paper [7] where the 
details of the proof have been written for several important cases, namely, Muskat, Water 
waves and SQG patches). 

Let us remark that, at the end, we have to work with a coupled system involving the 
evolution of the surface X, the "vorticity density" uj, the Ray leigh- Taylor condition a, the 
non-seffintersecting character of S quantified by the gauge F(X) and the tangential parts 
C\X ax + C2X a2 of the velocity field. 

This paper is a continuation of [5] where the two-dimensional case was considered. Many 
of the needed estimates can be obtained following exactly the same methods that where used 
in [5] for the lower dimensional case. Therefore, in order to simplify our presentation, we shall 
avoid here many details which were carefully proven in that quoted paper. This is specially 
the case of section 6 (control of the Birkhoff-Rott integral) , section 8 (energy estimates) and 
also for the approximation schemes which are identical to those developed in [5]. Therefore in 
the following we shall focus our attention on the more innovative parts of the proof, namely 
the evolution of the Rayleigh- Taylor condition, the non-selfintersecting property of the free 
boundary and the needed estimates for double layer potentials. 

5 Inverting the operator: The single and double layer poten- 
tials revisited 

Along this proof we need to consider the properties of single and double layer potentials, 
which are well-known characters in finding solutions to the Dirichlet and Neumann problems 
in domains D of W 1 . 

For our purposes those domains will be of three different types, namely: bounded, pe- 
riodic in the "horizontal" variables or asymptotically flat. We shall also assume that their 
boundaries are smooth enough (says C 2 ) and do not present self-intersections. Therefore one 
has tangent balls at every point of the boundary, one completely contained in D and the other 
in D c . We shall denote by v(x) the unit inner normal at the point x S dD, then under our hy- 
pothesis we have that, for r > small enough, the parallel surfaces dD r = {x+rv(x)\x £ dD} 
are also C 2 surfaces with curvatures controlled by those of dD. Furthermore the vector field 
v can be extended smoothly up to a collar neighborhood of dD allowing us to write the 
following formula: 

d tx du 

Au(x) = ^-^(x) - h(x) — (x) + A s u(x) 
ou z du 

where A denotes the ordinary laplacian in M. n , A s is the Laplace-Beltrami operator in dD, 
h(x) is the mean curvature of dD at the point x and u is any C 2 -function defined in a 
neighborhood of dD. 

For convenience we will use the notation D\ = D, D2 = D c , S = dDj and Vj(x)(j = 1,2) 
the inner normal at x 6 S pointing inside Dj. Let dS be the surface measure in S induced 
by Lebesgue measure in ambience space, then given integrable functions </?, ip on S we have 
the integrals 

V(x) = c n / ip(y) T , n77Z2 dS (y) 

Js If -y\\ 



W(x) = cj y(y)^-( .. 1 ~n^ )dS(y) 

Js ou x \\x — y\\ 

representing the single (respect, double) layer potential of ip (respect, ip), where c n is the 
normalizing constant so that , f^-'i becomes a fundamental solution of A in W 1 , n > 3. 

For x € S let us denote W\{x), V\{x) (resp. W2(x), T^(x)) the corresponding limits of the 
potentials inside D\ (resp. D%), we have: 

Wi(x) = -(<p(x) - J <p(y)K(x,y)da(y)) = -(<p(x) - V<p(x)) 
W 2 (x) = ±(<p(x) + J <p(y)K(x,y)da(y)) = ~{<p(x) + T><p(x)) 

f^(z) = -\(?l>{x) + J il>(y)K(y, x)da{y)) = -|(^(x) + V*^{x)) 

dV 1 f 1 

^-0) = -^WW - / i>(y)K(y,x)da(y)) = --(^(x) - D>(x)) 



where 



K{x,y) = 2cn-—{- 



' O VII I In— 9 / ft In 

ak'n |p — 2/|| F _ 2/1 

It is well known that in those scenarios considered above the boundary operators D (and 
P*)are smoothing of order —1 and therefore compact. Furthermore all their eigenvalues 
are real numbers having absolute value strictly less than 1. Therefore, by the standard 
Fredholm theory, the operators I — AT>, I — XD* are invertible when |A| < 1. However, in our 
case the domains are moving and the evolution of their common boundary S involves such 
inverse operators, making it necessary to estimate their norms in terms of the geometry and 
smoothness of S. 

Although there is a vast literature about single and double layer potentials, we have not 
been able to point out a precise statement giving the information needed for our results. 
Therefore in this section we provide arguments to prove that the norms of such inverse 
operators growth at most polynomially P(|||5|||), where |||S||| is just ||S||^2 plus a term of 
chord-arc type controlling the non-self-intersecting character of the boundary, namely we add 
the term r(5) _1 , where r(S) is the sup over all the positive r so that S admits tangent balls 



of radius r in both domains Df 



\\S\\ C2 + (r(S)Y 



We shall write P(|||5|||) to denote < C(|||5||| p ) for certain positive constants C,p which 
are independent of the characters whose evolution is being controlled, but the size of both 
constants may change along the proof and we shall make no effort to obtain their best values. 

We will consider the case of bounded domains in M n , n > 3, because the needed modifi- 
cations when n = 2, namely taking log |»| as fundamental solution for the laplacian, as well 
as the changes for the periodic or asymptotically flat domains, are left to the reader. 

Let T> and T>* be the potential defined above with kernel 

K( n _ _^ i _ {x-y,v{y)) 

J\ I JU • y I — ^n c\ / \ i \ ii o — ^72 i i 

ov(y) \\x — y||™ \x — y\ n 

and K{y,x) respectively. In the study of the inverse operators (I — XD)~ l , |A| < 1 it is 
convenient to consider first the particular values A = ±1. 
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Proposition 5.1 The following estimate holds 

\\(I±V)-i\\ L . {s) =P(\\\S\\\). 

Since the boundedness of (/ ± T>)~ 1 in L 2 (S) is well known from the general theory, 
we can simplify the proof considering only functions / G L 2 (S) whose support lies inside a 
region of S where the normal v(x) is close enough to a fixed direction. Then for a general 
/ an appropriate partition of unity would allow us to add the local estimates, so long as the 
number of pieces is controlled by |||/S|||. We shall use the following observation of Rellich 
(lemma I5.2[) whose proof is immediate 

Lemma 5.2 . Let u be a harmonic function and h a smooth vector field in the domain D, 
then we have: 

i) div (|Vu| 2 /i) = 2div ((Vn • h)Vu) + 0(|Vu| 2 |V/i|), 

it) / (v,h)\Vu\ 2 do: = 2 [ -^{Vu-h)da + 0{[ |Vu| 2 |V/i|). 
JdD JdD Of Jd 

Given a function / G C 1 (S) we may define V T / choosing at each point x € S an or- 
thonormal basis {ei, ...., e n _i} of the tangent space T X (S) (we can consider also V T f to be 
the gradient naturally associated to the induced Riemannian metric by the ambience space). 
In both ways, although different, we have that |V T /| = A T / is an elliptic pseudo-differential 
operator of order 1 in S. Solving the Dirichlet problem An = in D, u\s = / we obtain the 
operator D u = -^\s which is also a pseudo-differential operator of order 1 in S. 

Lemma 5.3 Let f G L 2 (S) having support on the region i < (v(x),rj) < 1 (for a fixed unit 
vector n), then we have: 

I \D u f\ 2 da ~ f \V T f\ 2 da 
Js Js 

where the constants involved in the stated equivalence ~ are P(\\\S\\\). 

Proof: Let u be harmonic in D so that u\s = f ■ Under our hypothesis about / and since 
|Vn| 2 = \D u u\ 2 + |V r n| 2 and V T u is a local operator (supps(V T /) C supp(/)), lemma 
yields: 

- / \V T f\ 2 da < / {u(x),i])\V T u\ 2 da < 3 / \D u u\ 2 da + 2 / \V T u\\D v u\da 
2 Js Js Js Js 

from which we easily obtain 



/|V r /| 2 ^<P(|||S|||) f \D u f\ 2 da. 
Js Js 



To get the opposite inequality we proceed as before, but since D u f is not local, an 
extra argument is needed to control the contribution of the region outside supp(/). Let us 
introduce surface discs B r (x) = {y G S\\\x — y\\ < r}, x G S, < r < |||5'||| _1 and domains 
A. r (x) = {y + pv(x)\y G B r (x),p < r}. Given R = 5|||5'||| _1 there exists a fixed unit vector 
n so that ^ — ( t/ (y), r l} — 1> f° r every y G Br(x) and also a smooth vector field h such that 
h\A R (x) = r), supp(/i) C A 2R (x) and %\h(x)\ < (h(x),v(x)), \\Vh\\ 2 < P( 

In order to obtain the estimate 



/|^/| 2 dt7<P(|||5|||) / \V T f\ 2 da 
Js Js 
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we may assume, without loss of generality, that supp (/) C -Br(x), for some x G S, and then 
prove that 

/ |A,/| 2 <icx<P(|||S|||) [ \V T f\ 2 da, 

JB R (y ) JS 

uniformly on j/o G £*• 

With the vector field /i defined above in A2fl(y) let us apply Rellich's estimate to get 

[ \D u f\ 2 (h,u{x))da(x) = [ (u,h)\V T f\ 2 da-2 [ D y fV T f ■ hda + 0( I \Vu\ 2 \Vh\) 
Js Js Js Jd 

where u satisfies Au = in D, u\s = f ■ We get easily 

! \D u f\ 2 (h, v{x))da{x) =0([ \V T f\ 2 da + I \Vu\ 2 \Vh\dx). 

JB R (y ) Js JD 

Then the proof will be finished if we can show that 



f |Vu| 2 |V/i|(ix < P(|||S|||) [ \V T f\ 2 da. 
Jd Js 



To see it let us consider the parallel surfaces S r = {x + rv(x)\x G S} (0 < r < \\\S\\\) and 
observe that 

u~(x + rv{x))da 



2 j / 2/ 

u da r ~ j 



and 



/ [u 2 (x + rv{x)) — u 2 (x)]da(x) = / / Vu 2 (x + tv{x)) ■ v{x)dtda 
Js Js Jo 

= 2 f u(y)Vu(y)-v(y)<2(f u 2 (y))l(f \Vu\ 2 (y))l 

where L r = {x + pv{x)\x G 5, < p < r}. 

Taking F{x + rv(x)) = f(x)X(x) {X= smooth cut-off) as a comparation function, Dirich- 
let's principle and Poincare's inequality give us the estimate 

/ |Vn| 2 < / |VF| 2 <C(f |V T /| 2 + / |/| 2 ) = 0([ \V T f\ 2 da). 
Jd Jd Js Js Js 

Therefore 

u 2 do r ~ I u 2 (x + rv{x))da < I f 2 (x)da + (f u 2 (y))^( i |V T /| 2 )i 

An integration in r, < r < R = IHS'IH -1 yields 

f u 2 dx<R[! f(x)da + (f u 2 {y))\{f \V T f\ 2 )^. 

That is 

f u 2 dx <CR I \V T f\ 2 da. 

J Lit J j 
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To conclude let us observe that 



f \Vu\ 2 \Vh\ = - [ Au 2 \Vh\ = - [ (Au 2 \Vh\ - u 2 A{\Vh\)) + - [ u 2 (\Vh\) 
Jd 2 J D 2 J D 2 J D 

= »/*!. | VA|4, -I //«*„+ 1/^V|/>| 
2J S dv 2J S dv 2J D 

<{f f 2 da)Hj\^\ 2 \Vh\ 2 da)h +C J f 2 da + C J u 2 . 
Proof of Proposition 15.11 As before let / E C 1 (S), supp (/) C Uq and let u be its single 



layer potential: 



u(x) = c n - --—rdS(y) 

Js \\x-y\r l 



Then taking derivatives on each domain Dj with respect to the normal direction and evalu- 
ating at 5 we get: 



^ = ~c/w -»•/(»)). 



By lemma [5^31 we know that 



dv l2 , f ._ l2 , f , dv l2 . 

\ da ~ / V-riM da ~ — — acr 



S dz'i 7 5 7 5 <9z/ 2 

where the constants involved in the equivalences ~ are all controlled by above (respect, 
below) by P(|||S|||) (respect. 1/P(|||5|||)). 

Since ^j- + ^j- = — / those estimates imply that 

1 



min(||/-P7||2,||/ + P7||2)> 



P(\\\S\\\) 



i.e. ||(/±P) _1 || = P(|||5|||). Then using an appropriate partition of unity, that estimate 
extends to a general / G L 2 {S). q.e.d. 

Next we shall consider Sobolev spaces H S (S), < s < 1, defined in the usual manner 
(i.e. throughout local coordinates charts). We have also elliptic pseudo-differential operator 
A s = (—A) 2 in such a way that 



I^ ( 5)^II/IIl 2 + I|a s /IIl 2 . 

Then H~ S (S) = (H S (S))* allows us to consider the negative case by duality under the 
pairing 

I (fnpda, <j> e H- s , ip e H s 
Js 

and 



H-s = sup / (pipda 
\M\ H s=iJs 

Since both T> and T>* are compact and smoothing operators of degree —1, the commutators 
[A S ,T>], [A S ,T>*] are then bounded in L 2 (S) (0 < s < 1) with norms controlled by |||5|||, 
allowing us to extend proposition 15.11 to the chain of Sobolev's spaces: 

Corollary 5.4 The norm of the operators (JzLP)" 1 , (ZiP*)" 1 in the space H S (S), —1 < 
s< 1, is bounded 6j/ P(|||5|||). 
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5.1 Estimates for (I + XV)- 1 , |A| < 1. 
With the same notation used before we have: 

1-XdV 1 + XdV \ tll . ^ lf .. 

1 + XdV 1-XdV l, ir . ^ mi/ ^ 

where 

V(x) = c n [ ^ dS(y). 

Js If- y\\ 



Then the identity 4> - XV* (p = yields 



= (1 - A) / V^-dS + (1 + A) / V^-dS = (1 - A) / | W| 2 + (1 + A) / | W| 2 

JdDx " v l JdD 2 dv 2 JDi Jd 2 

which implies eft = 0. Similarly for eft + AD*</> = 0, — 1 < A < 1. 

Remark: This observation can be improved applying the following fact (whose proof we 
skip because it will not be used in our theorem): 

|V«| 2 ~ / |Vu| 2 



D 1 JD 2 

where, again, the ~ is controlled by P(|||5|||). In particular it implies that the spectral radius 
of the operators V, V* is less than 1 - (P(|||S'|||))" 1 . 

Theorem 5.5 The operator norms \\(I + XV^Wh*^, ||(/ + AP*) _1 ||^ (S ), \s\ < 1, |A| < 1, 

are P(|||«S|||) (growth at most polynomially with \\\S\\\). 

Proof: The identity (I - V)~ l {I - XV) = I + (1 - X)(I - V)~ l V shows that the conclusion 
of the theorem follows easily when |1 — A| < pd i lg iii s and similarly when |1 + A| < p/ i ifgMh - 
Therefore, without loss of generality, we may assume that 

1-IAI> ' 



P(IIISIII)' 



Assume now that </> G H 2 {S) satisfies that \\<p\\ _i = 1 and 

\\<f>-XV*ct>\\ i < 



'B~* ~ P(\\\S\\\Y 
Then the single layer potential 

Js \\x-y\\ n 2 



satisfies the inequality 

,s V{ *- xv '« )dSl< -msw 

On the other hand one have 



/. 



I V(<j) - XV*<p)dS = (1 - A) I | W| 2 + (1 + A) I |W| 

Js Jdx Jd 2 
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implying the estimate 

// w + AI, ** )d ^ (1 + A) / Dl |w|2 + (1 - A) / D J w|2£ W' 

Then adding both inequalities together we would obtain 

l VHai mm 

which is impossible because of the following: 

Lemma 5.6 If V is the single layer potential of 4> then 

[ V{x)<t>{x)dS{x) = [ [ Hx) ^ n ] _ 2 dS(x)dS(y) > _^_|M|2 

Js JsJs\\x-y\\ P (\\\S\\\) H ?(S) 

Proof: Let us observe first that 

J f H^(y) da[x)da{y) = f _L |^(£)|^>o, 
JsJs \\x-y\\ n 2 7 K n |£r 

where 4>dS denotes the Fourier transform of the measure 4>dS supported on S. This implies 
that 

(^) = f f j^^dS( X )dS(y) 
JsJs \\x-y\\ n 



is an inner product satisfying: 



K^>|<<&$3<^>3 



and we wish to show that 

To see it let us observe first that given eft £ H~2(S) then its single layer potential u\S belong 
to the space Hi^S) satisfying: 



which can be proved easily using local coordinates. As a consequence we have 

<l>(x)<l>(y) jat-\iar..\ ^ r»nil olllMI jlii2 



JLL_dS(x)dS(y) < P(\\\S\\\)U\Y , 
! S J S \\x - y\\ n z H 2(5) 

In the opposite direction, since H~ s = (H s )* we have 

ll^lltf-* = SU P / 4>{x)f(x)do-(x). 
feH° Js 

Let us assume, for the moment, that given / £ H s there exists g G H s ~ l so that 

Js If -2/H 

and 1 1/| |jfs ±2 | |g| |#s-i (where we have used again the symbol m to denote equivalence modulo 
a factor P(|||5|||)). 
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Then 

\\</>\\h-° - SU P (4, 9), 
\\g\\ H s-i=i 

and taking 8=2,8 — 1 = — 5 we get 

||0|| H _i < P(\\\s\\\M<l>)*{9,9)* < i'dll^lllX^^'IMI^-j < p{\\\s\\\)(M)*- 

To close our argument it remains to solve the equation 

fW = C - / II g{y ln-2 dS (y) 

Js If -y\\ 

i.e. to prove that given given / G H s there exits g £ if s_1 satisfying the relation above. 
To see that let us consider the solution of the Dirichlet problem: 

( Au = in D 1 

\u\s = f 

and the equation 

-2^— = g-V g 

OV\ 

i.e. g = (I — V*)~ 1 (—2-0 L ). Then we claim that such g verifies the identity 



This is because the function 



II g{y ln-2 dS (y)- 

Is \\x-y\\ 



V ( X ) = C n [ || 9{y L-2 dS (y) 

Js \\ x -y\\ 



is harmonic in D\ and satisfies 



dV „* du 

-2—=g-V*g = -2— 



which implies that V = u in D\ and, therefore, taking limits up to the boundary we obtain 

/(*) = c - / n g{y l n - 2 dS(y). 

Js \\x-y\\ 

To finish the proof of theorem 15.51 let us consider for every r, < r < 1, the identity 

(/ - ap)" 1 a t = a t (/ - \vy l + (i - \vy l c T {i - xvy 1 

where the commutator C T = [DA T — A T T>] is a pseudodifferential operator of order r — 2 
whose bounds are controlled by IIISHJ. Then 



\\(I-\V)- l f\\ H s<\\(I-\Vr l f\\ H _ h +\\k^{I-XD)- l f\\ H _ h 

<h/ii i+iKi -A2?r i A-+i/|i 4 



< ll/H^ + ||A'+i/|| 1 <P(|||5|||)||/||^ 



q.e.d. 
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Remark 5.7 In the particular case of the sphere S = S n 1 (n > 2) the estimate of lemma 
5.6 becomes an identity: 



i ^^- 2 dS{ X )dS{y) = c n Uf 

gn — l J gn-1 | |X — y\\ tl ^ \p' 1 x ) 



; dS(x)dS(y) = c n \\4>\\ 

J gn—l J gn—1 | |X" — y\\'~ " 

for n > 3, and 



i 



log 1 1 a; - y\\(p(x)(j)(y)dS(x)dS(y) = c 2 \\4>\\ 



2 



1 



gljgl " r,, if-2(Sl) 

/or n = 2. 

Proof: We present the details when n > 3. TTie case n = 2 follows similarly. Let 4>(x) 
2~2 akYk(x) where Y k is a spherical harmonic of degree k normalized so that H^Hl 2 ^™- 1 ) — 
then we have 

t—* 2k + n — 2 H 2 (5) 

fc>i 



Claim: If k 7^ j then 



Y k (x)Yj(y) 

-^^dSixjdSiy) = 



is™- 1 Js™- 1 \\ x — y 

Taking the Fourier transform and using Plancherel we get 

Yk{x)Yj{y) Jof . Jaf . f 1 



is"- 1 Js™- 1 \\ x — y\ 

But it turns out that 



^±dS(x)dS(y) = f ^Y^dS(i)YjdS{0^ 

Jm. n I? I 



2 i?i 

where J u designs Bessel's function of order v, implying the claim. 
Therefore our estimate diagonalizes: 

— \Y k dS{i)\ 2 di = c / -| J k+ ^ (r)\ 2 dr 

l?r Jo r 



2 



°°W dr= 11 



and the following well known identity for Bessel's functions 

1 1 
2ji 

allows us to finish the proof. 

5.2 Estimates for f2 and u. 

In the following we shall consider asymptotically flat domains leaving to the reader the details 
of the periodic case. Since we have controlled the norms of the operator relating °. and X, 
we are in position to obtain the following inequality: 

Holla* < P(\\X\\l + \\F(X)\\l x + UIA^IM, (5-1) 

for k > 4, with P a polynomial function. Then Sobolev's embedding implies 

INI/** < P(\\X\\ 2 k+1 + \\F{X)\\l x + IHAr 1 !!^), (5-2) 
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for k > 3. We will present the proof (|5.ip when k = 4, because the cases k > 4 can be 
obtained with the same method. 
Theorem 15.51 in (|2.6p yields 

||0||^ = ||(/-A M P)- 1 (-2A p X 3 )|| H i<C7||(/-^P)^|| f/1 ||X3|| H i<P(|||5|||)||X3|| H1 

implying that 

||n|| H i < P(||X||1 + \\F(X)\\ 2 LX1 + IIIJVI^Hlcc). 

Next we will show that 

\\d 2 ai n\\ L2 < p(\\x\\l + \\F(x)f LOO + \\\N\- l \\ L ~)\\n\\ H1 (5.3) 

which together with the estimate for ||^||#i above will allows us to control d^ CI in terms of 
the free boundary. 

In order to do that we start with formula (|2.8|) to get d? CI = 1^ + 1 2 + 13 + 14 — 2A p d 2 (1 X 3 
where 

2vr 7 R2 |X(a) -X(a-/3)| d 

2vr J R 2 \X(a) - X[a - p)\ A 

4vr J R 2 \X[a) - X(a- p)\° 

with A(a,/9) = (X(a) - X(a - £)) • (d Ql X(a<) - d ai A~(a - £)), and 

2vr J R 2 |A(a) - A(a-/^)|' :i 



Our next objective is to introduce the operators Tk (|9.5p defined at the appendix, in the 
analysis of those integrals Ij. Formula f|2.3[) gives us to = d a2 (Cld ai X) — d ai (Cld a2 X) and 
from standard Sobolev's estimates we get 

\\ij\lv < p(\\x\\l + ilF^nicc + |||jvr 1 || L ao)||n||j f i, j = 1,2, 

and similarly with I3. 
Regarding 

h= [ dp+ I dp = J 1 + J 2 
J\p\>i J\P\<1 

we integrate by parts in J\ to obtain 



I/3|>1 

A„ /• X(a) - A~(a - /3) „, „. n , n TW . 

27r 7^1=! |A (a) - A (a - p)\ 6 

From this last expression it is easy to deduce the inequality 

J 1 <C||F(X)||| oo ||A--(a,0)|| 2 cl ( / ^ a -/^ dp+ I \u(a-P)\dl{P)) 

J\p\>i \P\ J\p\=\ 

providing us with an appropriated control (see appendix for more details). 
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Next let us consider J 2 = K\ + K 2 + K% + K^ where 

Ki = ^PV [ >*f a )~*l a ~f!L A d a Ma ~ P)d 2 ai X(a - p)dfi • d ai X(a), 
27r J\p\<i\ x {<x)- X {.<x-P)\ 

K2 = ^PV [ |5 (a ?~w a ~f?i3 A d ai d a2 Q(a ~ P)d ai X(a - fid? ■ d ai X(a), 

27T J\/3\<i\X{a) - X{a- p)^ 

K 3 = ~^PV [ 3 a }~Y( a ~S\3 A a «i n ( Q - P)d ai d Q2 X(a - /3)d/3 • d ai X(a), 
27r J\p\<i \X{a,)-X(a- p)\ 6 

27T J\p\<i \ x { a ) - X \a-P)\ 

Then the terms K\ and Kj, are handled with the same approach used for I 2 ( i.e. f|9. 13|) in 
the appendix) and we rewrite K 2 in the form 

K 2 = ^f 3 a) ,~^ a ~oL A d ai d a2 £l(a-P)(d ai X(a-l3)-d ai X(a))dl3 ■ d ai X(a), 

to show that it can be estimated via an integration by parts in the variable (3\ using the 
identity 

d ai d a2 n(a- (3) = -d^{d aa ^l{a-0)) 

and the fact that the kernel in the integral K 2 has degree —1. 

It remains to deal with K4: to do that let us consider K^ = L\ + L 2 where 

Li = ^PV[ ,y/ a ?~y/ a ~/f?i3 A %M<*-P)( d °* x (.<x)- d c* x (<*-0)W ■ 9 ai x (a), 

27r J\p\<i |A(a)-A(a-/^)| d 

and 

2vr J\p\<\\ x \<x)- X k<x- P)\ 
The term L\ can be controlled like K 2 , and L 2 can be rewritten in the form: 

L2 = ^^ i,<i l |X(a)-X(a-/3)|3 " WxJaTWr aMa " W ' ^ 



showing that it can be estimated as we did with Ta (I9.8p . that is we obtain f|5.3j) . Similarly, 
equation (|2.9p yields 

||^ 2 n|| i2 < p(\\x\\l + \\F(x)f LOO + \\\N\- l \\ L ~m\\ m , 

and then the inequality 2||(9 Ql (9 Q2 ri|| L 2 < ||^ ftH^a + || d^ 2 Q\\ L 2 gives us the desired control 
upon [|n[|#2. 

Next we will show that 



^|| L2 < p(||x||1 + \\F{x)\\ 2 LOO + niivr 1 !!^)^!!^ (5.4) 



allowing us to use the estimates for ||f&||#2 above. In order to do that we start with formula 
(|2.8p to get d% ll Q = d ai h + d ai h + d ai l3 + d ai I4 — 2 A p d^ 1 X3 where the most singular terms 
are given by 

J * = ¥ pv \ 1 x(1 " x( a "ml A *<* ~ W ■ d ^ x ^ 

2vr J R 2 \X(a) -X(a- fj)\ A 
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A^ py r d 2 ai X(a)-d 2 ai X (a-f3) 
2ir J R 2 \X(a) - X(a 



^ = ^PV I . a \ v L, VFZ «M3 A w (« - W ■ d ^ X ^)^ 



3 ^M m , f t,,- m X («) " X (« " P) 



h = -V'PV / B(a, P) ' -^ A w(a - /3)d/3 • ai X(a), 

withfi(a,/3) = (X(a)-X(a-/3))-(a2 i X(a)-a2 i X(a-/3)), and 

2vr J R 2 \X(a) - X{a-p)\ 6 

and where the remainder terms can be estimated with the same method used before. 
Now we write 

^3 = -^Ti{d a2 (tid ai x) - d ai (ttd a2 x)) • d 3 ai x 

to obtain: 

\\j 3 \\ L 2<c\\Ti(d a2 (nd ai x)-d ai (nd a2 x))\\ L 4^ ai x\\ L4 . 

Next let us observe that in the proof of estimate (|9.9p one can replace L 2 by L p for 1 < p < oo 
(see |17j). In particular we have 

||J 3 || L 2<P(||X-(a,0)|| cl ,*+||F(A)|| L c*,+|^ 

and then Sobolev's embedding in dimension two: (||<7||jr,4 < C||g||#i) yields the desired control. 
Regarding J4 we follow the approach taken before for 73 but using now the L 4 norm. That 
is we split 



J A = / dp + / d(3 = K 5 + K 6 

J\fS\>l ^|/3|<1 

and since 

\u(a-P)\ 



K 5 < \\X - (a,0)\\c4 F (X)WU / ' \,|3 — W 

J\p\>i \P\ 

that term can be estimated as above. 

Next we introduce the splitting Kq = L3 + L4 where 

Ls = ^ i\d 2 ai X(a)-d 2 ai X(a-P))[ \ * ]Au(a-P)dp-d ai X(a), 

2nJ\p\<i \X{u)-X(a-P)\ A \VX(a)-P\ 6 

K f d 2 X(a)-d 2 X(a-P) , 

La = ^PV / ai \>. f 3 P) A u(a - P)dp ■ d ai X(a). 

2vr J m<l \VX(a)-P\ 6 

We have 

L 3 <C\\X- (a,0)\\M\\nX)\\i- +\\X- (aMlUWr'WU) f '^I'fV 



/9|<1 



W 



(see appendix for more details), giving us the appropriated estimate. Regarding L4 we use 
identity f|9. 16|) which after a careful integration by parts yields 

K f \P\(d 2 ai X(a) - d 2 ai X(a - p)) A uja - P) ■ d ai X(a) 
~27rJ m=1 |VX(a)-/3|3 d/ ^ 
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helping us to prove the inequality 

\\L4 L 2 < P(\\X - (a,0)\\ C 2 + \\F(X)\\ LOO + |||iV|- 1 || i o )(||^ 1 X|| L 4||a;|U4 + || W || L2 ). 

Clearly J5 can be approached with the same method used for J4. Regarding the term Jg we 
have to decompose further: first its most singular terms which are given by 

27r J\/3\<1 \ x { a ) -X(a-p)\ 6 

L 6 = ^PV f ly/^"^""! A 8 2 ai d a2 n(» ~ P)d ai X(a - p)d? ■ d ai X(a), 
27T J\p\<i \X[a) -X{a-p)\ a 

L 7 = -^L PV f X(a)-X(a-^) a _ 2 _ Q 

27r J\0\<i |A(a)-A(a-/3)| d 

L 8 = -^PV I ,y/ a ?~y/ a ~f?| 3 A %M<* ~ P)d* 2 X{a - p)dp ■ d ai X(a). 
27r J\/3\<i \X{a) - X(a-P)\ 6 

Second let us observe that the remainder is easy to handle; the terms L5 and Lj can be 
estimated as we did with K\ and K 3 using the L 4 norm and, finally, Lq and L$ are like K2 
and K4 respectively. Putting together all these facts we obtain (|5.4p . 
Similarly to the case of lower derivatives, equation (|2.9p yields 



\\n\\ H3 < p(\\x\\l + \\f(x)\\ 2 loo + niivr 1 !!^)^!!^. 

To finish it remains to show the corresponding inequality for derivatives of fourth order: 

\\n\\ m < p(\\x\\l + UFpoiiico + iWN^WL^wnWHS. (5.5) 

Identity (\2.8h allows us to point out the most singular terms in d^ fl: 

M ' = ^ PV f 1 xm " xt "ml A ^ a ~ PW ■ W a )» 

2vr y H 2 \X(a) - X(a-p)\ A 



M = A Epv [ dl^-dlXja-p) 
2 2tt y M2 \X(a)-X(a-p)\* 



A oj(a - p~)dp~ ■ d ai X(a), 



M 3 = - 3 -^PV [ C{a,P) ^ a) * (a ^ 5 Au(a-{3)dp--d ai X(a) 

4vr J R2 \X(a) - X{a - j3)\ b 

with C{a,{3) = (X{a) - X(a - (3)) ■ (d^Xia) - d^Xia - p)), and 



2vr J R 2 \X(a) - X(a 
Then in order to estimate Mi we start with [|Mi||^2 < Cif||9 4 X||^2 where 



K = sup 

a 
Following ref. [8] we have: 



™l vtS-Z-w **- 1 ™ 



K<0 1 + 2 + 3 + 4 + 5 
where 



0\ = sup 

a 
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2 = sup 



X(a) - X(a - (3) - VX(a) ■ /3 



O3 = sup 



/9|<1 
VX{a)-/3{- 



\X{a)-X(a-f3)\ 3 
1 



Au)(a- f3)dp 



1 



l/3|<l 
O4 = sup 



X(o)-X(a-/3)|3 |VX(a)-/3|3 
VI(a) • /3 



A oj(a - j3)d/3 



5 = sup PV 

a 

An integration by parts in 0\ yields 

O x <C , ||VX||| 00 ||F(X)||| 00 sup( 

a 

<C\\VX\\ 2 Laa \\F{X)f L oo\MLo 



m<1 \VX{a)-P\* 

VI(a) • p 



p\ <x \VX(a)-p\ 

Ma-(3)\ 



A (w(a- /3) -w(a))d/3 

„ A u(a)d[3 



I/3|>1 



|/3| 



d/3 + / |fi(a-/3)|dZ(/3)) 

Vi=i 



and Sobolev's embedding allows us to conclude. 
Regarding O2 we have 

2 < \\X- (a, 0)11^11^(^)111- llwlUcc I / |/3| 2 ~ 5 (i/3| 

J|/3|<1 

and then the estimate, ||w|| c « < C||u;||#2 for < 5 < 1, gives the desired control. Using 
(I9.15P and after some straightforward algebraic manipulations we get a similar inequality for 
O3. Next we have 



O4<C||X-(a,0)||^ 1 |||iV|- 1 ||| O o||w|| c ,| 



|2-5 



dP\, 



I£I<1 



giving us also the same estimate. Furthermore it is easy to prove that O5 = 0. 
Next we consider the term M2 with the splitting: M2 = Qi + Q2 + Q3 where 



Oi- 



d 3 ai X( a )-d 3 ai X(a-(5) 
2vr ./|/3|>i \X{a) - X(a - P)\ A 



Q2 = ^ 



d 3 ai X(a)-dlX(a-(3) 
2xJ m<1 \X(a)-X(a-P)\ 3 



A (a; (a - /?) - oj(a))df3 ■ d ai X(a), 



d 3 ai X(a)-d 3 ai X(a-P) 



= to PV J w<1 \X(a)-X(a-[3)\ 3 



dj3 A a; (a) • <9 Ql X(a). 



The term Qi can be estimated as before, regarding Q2 we can use the identity 



d^Xia) - d^X(a - /3) = / Vd A at X{a + (s - l)/3)ds ■ f3 

Jo 



and the control of Q3 can be approached as we did with the operator in (J9.7)) . Similarly with 
M3, whether M4 is analogous to Jq, and all these observations together allow us to obtain 
(ESD. 
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6 Controlling the Birkhoff-Rott integral 

Here we consider estimates for the Birkhoff-Rott integral along a non-selfintersecting surface. 
Let us assume that V(X(q) - (a,0)) G H k (R 2 ) for k > 3, and that both F(X) and jiVl" 1 
are in L°° where 

F{X){a,P) = \P\/\X{a)-X(a-P)\ and N{a) = d ai X(a) A d a2 X{a). 

The main purpose of this section is to prove the following estimate: 

\\BlUX,ui)\\ H *-i < P(\\X\\l + \\F{X)f Loa + IHAr 1 !!^), (6-1) 

for k > 4. Here we shall show it when k = 4, because the other cases, k > 4, follow by similar 
arguments. We rewrite BR in the following manner: 

BR(X,u)(a,t) = ~FV f * (a) ~ *^ A (dp 2 (nd Pl X) - d^(Qd^X))(P)dp, 

4vr J M 2 |A(aj - X(P)\ 6 

which together with the estimates about f2 in section 5 and also about the operator 71 in the 
appendix, yields 

\\BR(X,u)\\ L2 < P(\\X\\l + \\F(X)f Loa + IHiVr 1 !!^). 

To estimate derivatives of order 3 we consider d^.(BR(X,uj)), and observe that the most 
dangerous terms are given by 

T J_ pT/ f KX{a) - d^Xja - 0)) A u(a - g) 

11 ~ A* L \X(a)-X(a-P)\* ^' 

I 2 = —PV (X(a)-X{a-P))Auj(a-P) ' — -^ dp 

4vr / E 2 |X(a)-X(a-/3)| 5 

. J_ py /" {X{a) - X{a- P)) /xjdl^ja- P) 

3 4tt y M2 |A-(a)-X(a-/3)|3 P ' 

In the appendix we find all the ingredients needed to estimate these terms Ij while the 
remainder in d^.(BR(X,u)) is easily bounded, namely: in 1$ we can recognize an operator 
with the form of 71 in (J9.5J) . so the estimate for uj in section 5 gives the desired control for 
I3. Regarding I\ we may use the splitting I\ = J\ + J 2 where 

J_ /• (^ A-(q) - dlX(a -P))A (q,(q) - q,(q - /?)) 
1 4vri E2 |X(a)-X(a-/3)|3 "' 

J2 -^T APV L \X(a)-X(a-P)\* ^ 

Then the identity fig. X(a) - fig. X(a - P) = p • J* Vd^X(a + (s - l)p)ds allows us to find in 
Ji a kernel of degree —1 which we know how to handle (see appendix). One use the estimate 
for TJj (|9.7p to deal with J2 and we proceed similarly to control 1%. 
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7 In search of the Rayleigh-Taylor condition 

As it was pointed out in section 4 (outline of the proof) our approach is based on energy 
estimates and a crucial step is to characterize those terms involving higher derivatives which 
are controlled because they have the appropriated sign. In our terminology they constitute 
the Rayleigh-Taylor condition, which is supposed to holds at time T = 0, being an important 
part of the proof to show that it prevails under the evolution. 
Let us introduce the notation 

\\\x\\\ 2 = \\x\\ 2 + \\F(x)\\ 2 LOO + 111x1-%^ 

where 

||X|| fe = || Xi - aillxa + \\X 2 - a 2 || L3 + ||X 3 || L2 + \\V{X - (a, O)))!^ , (7.1) 

and 
||V(X - (a,0))ll2r*-i = l|V(X - (a,0))||| 2 + \\d k ai (X - (a,0))||| 2 + \\d k a2 (X - (a,0))||| 2 . 
In order to justify the formula 

+ P{\\\x\ut)\ 

(here k > 4, although for the sake of simplicity we will present the explicit computations when 
k = 4, leaving the other cases as an exercise for the interested reader), it will be convenient to 
make use of the following tools, giving us different kind of cancelations, and which constitute 
our particular bestiary of formulas for this paper: 

From the definition of the isothermal parameterization we have the identities: 

\8 ai X\ 2 = \d a2 X\\ (7.2) 

d ai X-d a2 X = 0, (7.3) 

which yield 

^A(\d ai X\ 2 ) = \d ai d a2 X\ 2 - d 2 ai X ■ d 2 a2 X, (7.4) 



d^x ■ d ai x = -3d 3 ai x • d 2 ai x + (dl^- l d ai ){\d ai d a2 x\ 2 - d 2 ai x • d 2 a2 x), (7.5) 
di 2 x ■ d a2 x = -3d 3 a2 x • d 2 a2 x + (d 2 a2 ^ l d a2 ){\d ai d a2 x\ 2 - d 2 ai x • d 2 a2 x). (7.6) 

Using (I7.3J) and (J7.4J) we obtain: 

)r) 3 Y . r) f) Y - f) 2 f) Y . Ft 2 Y 

(7.7) 



d 4 ai x ■ d a2 x = -2d 3 ai x ■ d ai d a2 x - d 2 ai d a2 x ■ d 2 ai x 

- (d ai d a2 A~ l d ai )(\d ai d a2 x\ 2 - d 2 ai x • d 2 a2 x), 



di 2 x ■ d ai x = -2dl 2 x ■ d ai d a2 x - d 2 a2 d ai x • d 2 a2 x 



(d^d^A-'d^dd^d^x] 2 - d 2 ai x ■ d 2 2 x). 



(7.8) 
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And Sobolev inequalities imply that if V(X - (a,0)) G # 3 then d^X • 3 aj A~ G F 3 for 
i,i = l,2. 

With the help of the estimates above we may now afford the task of determining a. There 
is a part that may be considered as a mere "algebraic" manipulation to detect the relevant 
characters and, in so doing, we disregard many terms because they are of lower order in the 
sense of Sobolev spaces. At the end, we shall present how to deal with those lower order 
terms, if not for the whole collection of them, at least for the ones that we may consider 
to be the most "dangerous" characters. Here it is convenient to recommend the reader our 
previous works [H [5] where similar estimates were carried out. 

7.1 Low order norms 

Since -Xj(a) — > a.% for i = 1, 2 at infinity, let us consider the evolution of the L 3 norm. That 

is 



- — \ 



where 



Xi-ai||| 3 (i)= / \X 1 -a 1 \(Xi-a 1 )X lt da = I 1 +I 2 + h, 

h = \X\ — ai\(X\ — oti)BRida, 

h = / |-^i — ai|(^i — ai)Cid ai X\da, h = \X\ — ai\(X\ — ai)C2d a2 X\da. 

Then we have 

h < \\X X - ai \\l :i \\BR\\ L3 < C(\\Xi - axilla + H^IMI^H^), 

and Sobolev estimates, together with (|6.1|) , yield the appropriate control in terms of P( \ \ \ X \ \ \ & ) . 
Next since d ai X\ — > 1 as a — > oo, we have 

h < ||9o:i-X"l||z,°° ||-X"l _ "1 llz,3 || C*! ||£,3 , 

and it remains to get control of C\. Using (|3.ip we introduce the splitting C\ = Yl l j=i^i-i 
where 

2ir J R 2 \a-(3\ 2 \Xp 2 \ 2 2vr J R2 \a - (3\ 2 \Xp 2 \ 2 



2tt J R2 \a - f3\ 2 IXfcl 2 2vr J R2 \a - /3\ 2 \X^\ 2 

We shall show how control C\, because the estimates for the other terms follow by similar 
arguments. Integrating by parts one obtain C\ = D\ + D 2 where 



2vr J R2 \a-P\ 2 H \\Xfo\ z J vr J R2 |a-/3| 4 \Xp 2 \ 2 

Regarding D\ we write D\ = E\ + E 2 where 

E i = T~ [ ^BR(a-P)-dJ-§^)(a-P)d(3, 

27T J\/3\<1 \P\ y \ X P2\ J 
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27r J\3\>1 \P\ K \ X 2 \ ' 



2vr^| >1 |/3| 2 v "' M2 V|X/3 2 | 
Then Minkowski and Young inequalities yield respectively 



lEkU* < C\\BR-dJj^)\\ LS < P(\\\X\\U) 



\\E 2 \\ L , < C\\BR.dJ-§^)\\ L1 < C\\BR\\ L2 \\dJ-^)\\ L , < P(\\\X\\U), 

and the desired control is achieved. In the term D 2 we have a double Riesz transform and 
the standard C alder on- Zygmund theory yields 

\\D 2 \\ L z < C\\BR-^%\\ L3 < C\\\X^\- X \\ L ^\\BR\\ L3 < P(|||X||| 4 ). 
\ X M 

The estimate for ^3 follows on similar path, and the case of the second coordinate is also 
identical: 

^||^ 2 -a 2 || 3 L 3(t)<P(|||X||| 4 ). 

Regarding the third coordinate we have a stronger decay because of the asymptotic flatness 
hypothesis: 

li\\X3\\h(t) = f X 3 BR 3 da + [ X 3 dd ai X 3 da+ [ X 3 C 2 d a2 X 3 da 

* at JR2 Jr2 J K 2 

X 3 BR 3 da -\ I (d ai d +d a2 C 2 )\X 3 \ 2 da, 

2 2 JR2 



therefore the use of Sobolev's embedding in the formulas for C\ (|3.ip and C 2 (|3.2p . together 
with the estimates for BR (|6.ip . allows us to obtain: 

\j t \\x 3 \\Ut)<p(\\\x\\U). 

Once we have control of higher order derivatives, we can use the estimates of the appendix 
to get 

iA||v(X-(a ) 0))||| 2 (i)<P(|||X||| 4 ). 

7.2 Higher order norms 

Let us consider now 

\j t W^X\\h(t) = Jj^X ■ d i ai BR(X,co)da 

+ [ d A ai x ■ d A ai {dd ai x)da + f d A ai x ■ d A ai (c 2 d a2 x)da 

JR 2 JR 2 

= h + I 2 + I 3 , 

The higher order terms in I 2 and ^3 are given by 



(7.9) 



J\ = I C x d A ai X ■ d 5 ai Xda, J 2 = I 8LX ■ d^Xd^dda 

Jr 2 Jr 2 
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h = I C 2 d A ai X ■ d A ai 3 a2 Xda, J 4 = i d A ai X- d a2 Xd A ai C 2 da 

JR 2 JR 2 

Integration by parts yields 

Jx + h = -\ [ (d ai Ci + d a2 C 2 )\d A ai X\ 2 da 
2 Jr 2 

and therefore 

J\ + h < ^(l|5 ai Ci|| Loo + 11^(7211^)11^X11^ < P(|||X||| 4 ). 
Then in J 2 we use (J7.5J) to get 

h = - [ d ai {d A ai x ■ d ai x)dl 1 c l da < \\d ai (d A ai x ■ d^WvWd^dWv. 

Jr. 2 
Whether in J4 we use (|7.7|) to obtain 

Ja = - I d ai {d A ai x ■ d a2 x)dl 1 c 2 da < \\d ai (d A ai x ■ d a2 x)\\ L 2\\d 3 ai c 2 \\ L 2. 

Jr 2 



From formulas (|3.ip , (|3.2p one realizes that C\ and C 2 are at the same level than Birkhoff-Rott 
(|2.5|) . and, therefore, we can use the estimates for BR (|6.1|) to control ||flg C* 1 1 z, 2 » i = 1,2. 
Then formulas (|7.5p and (|7.7p indicate how to estimate \\d ai {d A il X ■ d ai X)\\ L 2, i = 1,2. That 
is we have: 

■h + JA<P{\\\x\\U). 

In ii the most singular terms are given by 

47T ./ib2./tb2 Alal-AltfiP 



3 /" /" 4 (X(a)-X(/3))-(9i X(a)-0* A7/3)) 

47T ./iB2./ro2 Ala)— JL D)P 



(7.10) 
1 pt/ /" /" * y, , (A'(a)-X(/3))A^ 1 a;)(/3) j 

47T ./ib2 ./iD.2 A (a) - A « 



Let us consider now the splitting J5 = K\ + if 2 

If! = -i-FK / / d A ai X(a) A (^X(a) - ^ X(/3)) • " (/?) + *ff dad/3, 
wr Jir2 7 R 2 |A(a) - A(/3jp 

K 2 = ±PV [ [ d 4 ai X(a) A (^X(a) - B^Xtf)) • " (a) ~ y^L ^A 
ott Jr2J R 2 |A(a)-A(^)p 

Next we exchange a and /3 in K\ to get 

Ki - h pv L L Bi ' xw) A <8 - x(a) - 8 ^ xm ■ \xlT-xw daiP 



167T 
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and therefore we can conclude that K\ = 0. In K 2 we find a singular integral with a kernel 
of degree —2 

K 2 = -±PV [ S^X(a) ■ f d 4 ai X(P) A "W~"W dflda, 
ott J R 2 J R 2 \X(a)-X{P)\ a 

and as it is proved in the appendix we have 

K 2 <P(\\\X\\U). 
Let us now decompose Jq = K% + K\ + K\ + K\ + K§ where 

K 3 = —PV / / d^X( a ) ■ (X(a) - X(0)) A u[fi) ——\ ,, |5 1 dad/3 

4vr J R 2 J R 2 \X(a) - X{fi)\° 

with A(a, (3) = X(a) - X(/3) - VX(a)(a - /3), 

TCi " 3 PT/ f f M Y( \(X< \ TT(R\\K { R\ ^ Z fr^gX (a) - d a X (J3)) • d^Xj/3) 

K 4 = —PV / / d ai X(a)iX(a)-X{P))Au{P) — — — — 1 dad/3 

** 3 pv /" /V w w*7 ^ r/fflu .m MjKW) • <X(a)-^X(/3) • <X(/3)) 

^ 5 = 4^ L lSX(a)iX(a)-X{P))^(l3) \X(a\ - X(»™ "^^ 



In K% and 7Q we find kernels of degree —2 and, as it is shown in the appendix, they behave as 
a Riesz transform acting on d^ X. In K\ the kernels have degree —3 and act as a A operator 
on d a X • d 4 il X. Then using the formulas (17. 5j) and (17. 7j) we get finally the desired estimate. 
We will find the R-T condition in J7. Let us take Jj = Kq + K7 where 

1 D¥ / „ 4 yn ( ,(X(a)-X(P)) VA»(a-/3) , _ 4 um ,„, 

^6 = -;r PF / d ai X{a) ■ / 1 V/ / y/mia " i uy / u mW) A {^){P)dpda, 

4vr 7 R 2 J R2 v |X(a) -X(/3)| d |VX(a)(a - /3)| J/ 

K 7 = ~PV f d 4 ai X(a)- f Jyff/" " f?| 8 A {fa U ){p)dpda. 

4vr y R2 y K 2 |VA(a)(a-^)| d 

The term Kq is controlled by (|9.8p in the appendix. Using (|7.2p and (|7.3p we get 

^7 = ~~PV lrl ai J a L • ^ J («) A *i(^i<")(«) + ^ X («) A R2(di 1 u;)(a))da. 
2 J R 2 \d ai X(a)\ 6 

Formula (|2.3p help us to detect the most singular terms inside K7, which will be denoted by 
Li, i = 1, ..., 8 and are the following: 

L x = -\PV f d 4 ai X(a) ■ ® ai ^% A R l (d 4 ai d a2 nd ai X)(a)da, 

L 2 = ~\PV I d^X{fx) ■ ,.? ai y/-l A R^nd^X^da, 

L 3 = ipy / ^X(a) • f ai ^ (a) A R 1 (d 5 ai nd a2 X)(a)da, 
2 Jr2 |d ai A(a)| d 

L 4 = \pV I 6f^X{a) ■ |f ai w a l /\Ri(d ai nd 4 ai d a2 X)(a)da, 
2 JR2 |CaiA(aJ|° 

U = ~\PV I d^Xia) ■ ^ A ^(^^na^XJCaJda, 
2 7r2 |a a2 X(a)| d 
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U = ~\PV I S^X{a) • i^yffia A R 2 (d a2 nd 5 ai X)(a)da, 

L 7 = ipy f d^Xia) • i^y/"^ A P 2 (^d a2 X)(a)cfo, 

L 8 = ipy / ^X(a) • MiAfi 2 (9 ai <9 tt2 I)( tt )d Q . 
In L\ we get a kernel of degree —1 of the form 

Ll = l -PV [ &^X(a) ■ [ r "fa |f ttl y/ a ?| 3 A (9 ai X(a) - ^X^^O^cta, 

2 Jr2 J r2 |a-/3| d |5 q , 1 X(q ! )|' j 

which can be estimated integrating by parts throughout d^d^O,; also the term L-j follows in 
a similar manner. In order to estimate L 2 , L±, Lq and Lg we realize that they can be written 
like (j9.3|) in the appendix plus commutators of the form (|9.1|) . Next we have to deal with L 3 
and L5: With L 3 we proceed as follows 

l 3 <l 3 + iWd^x^w^Wd^xw^WR^nd^x) - R^tydonXWv 

where L 3 is given by 

L 3 = \PV J d 4 ai X(a) ■ N ^ ] (R 1 d ai M i n)(a)da, (7.11) 

2 Jr2 \d ai X[a)\ 6 



and the commutator estimates (|9.ip show that it only remains to control L 3 . We use now 
formula (j!T8j) to get Z 3 = M x + M 2 where 

Mi = -A p PV f d 4 ai X(a) ■ N W (R 1 d ai )(d 4 ai X 3 )(a)da, 

J R 2 |d Ql X(a)| d 



and 



M 2 = -A^PV [ d^Xia) ■ N ^y (Rxd^ ) (^ (£fl(X, u) ■ d ai X))(a)da. 

JR2 \o ai A(a)\° 



Then we write Mi = 0\ + 2 + O3 where 

91 X 



o x = -a p pv / -^_L (a ai x 2 a a2 x 3 - d ai x 3 d a2 x 2 )(R 1 d ai )(dt 1 x 3 )da, 

JR 2 |C Q iA| 

2 = -A p Py / |/ 1 y| 3 (^ 1 X 3 a2 Xi - a ai Xi5 Q2 X 3 )( J Ri9 Ql )(^ 1 X 3 )da, 

Jr2 |a ai A|° 

3 = -A p PV i .. N3 v]3 di 1 X 3 (R 1 d ai )(di 1 X 3 )da. (7.12) 

JR 2 \0 ai Ji\ o 



Next we consider 0\ = P\ + P 2 + P 3 with 



P 1 = -A p PV [ . " 1 ^ 9 ai X 2 (Pi9 ai )(a a2 X 3 9^ 1 X3)(ia, 

/• a4 y- 

p 2 = a p pf / ai I a Q2 x 2 ( J Ria ai )(a ai x 3 ^ 1 x 3 )d«, 

JR2 \O ai JL\° 
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P 3 =A p PV f j£^d ai X 2 [(Rid ai )(d aa X 3 d^X 3 ) - d^X^R.d^id^X^da 

j R 2 |a Ql A| 

+ A p PV / T^ L ^d a2 X 2 [d ctl X z {R l d ai ){di 1 X z ) - {R x a ai ){d^X^X^)]da 

and the commutator estimate allows us to control the term P3. 
Now we use (jT.Tj) to write P\ = Q\ + Q 2 + Q3 

Q 1 = A p PV [ .J* 1 ] d ai X 2 (R 1 d ai )(d a2 X 1 di 1 X 1 )da, 

Jr2 \o ai Ji\° 

Q 2 = A p PV [ ai *] d ai X 2 {R 1 d ai ) (8 a2 X 2 8^ X 2 )da, 

Jk 2 \d ai ^\ 

Q 3 = A p PV / ai „| 3 ai X 2 CRia ai )(l.o.t.)cfa. 

JR2 |C Ql A|° 

The term Q3 is easily estimated. Regarding P 2 equality (j7.5[) allows us to write P 2 = Q\ + 
Q5 + $6 where 

Q 4 = -A p Py / ai ] d^XziRidc^id^Xi^X^da, 

Jm 2 |C Ql A|° 

r 8^ x 
Q 5 = -ApPV / . a ai v ], d a2 X 2 (R 1 d ai )(d ai X 2 di 1 X 2 )da, 

JR2 l^axAl 

/" 8^ X 
Q 6 = -ApPV / .J" 1 Y 1 n d aa X 2 (Rid ai )(lo.t.)da. 

Let us recall the identity Pi + P 2 = (Q 4 + Qi) + (Q 2 + Q5) + (Q3 + Qg) where (J 3 and Q6 
are easily estimated. With respect to Q 2 + Q§ we have 

Q 2 + Q b = ApPV f J ^^-d ai X 2 [(R l d ai )(d a2 X 2 di 1 X 2 ) - d a2 X 2 (R 1 d ai )(8i 1 X 2 )]da 
Jm 2 \d ai ^\ 

+ A p PV I jTr^d^X^d^X^Rid^id^X^ - (R x d aj )(d ai X 2 a^ 1 X 2 )]da 

Jm 2 |Cai A r 

and again the commutator estimates yields the desired control. 
Next we have 

Q 4 + Q x = ApPV f Ql ] d^X^X^R^id^Xi) - (RidaJid^X^Xifida 

JM. 2 I^ai^-I 

f 8^ X 
+ A p PV / i a a 'J, 9a 1 X2[(fiia ai )(a a ^iCli) - 5 Q2 Xi(i?i5 Q1 ) (5^Xi)] da 

Jw 2 |C Q1 A| 

~ A p PV [ T7^ F T 3 d i ai X 1 (R 1 8 ai )(d 4 ai X 1 )da. 

Jm 2 \o ai Jt\ 

The first two integrals above are easily handled allowing us to get 

Oi = P l +P 2 + P 3 <P(\\\X\\\ A )-A p PV I -^—di 1 X 1 (R 1 d ai ){dt 1 X 1 )da. (7.13) 

Jm 2 \o ai X\* 
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For the term 2 we proceed in a similar manner, first we check that 2 = Pa + P5 + Pq 

f ffl X 

P 4 = A p PV / ,^ ai 2 d ai x l (R l a ai )(d a2 x 3 di 1 x 3 )da, 

Jm 2 |Cai A r 
P 5 = -A p PV I rf^d^XxiRtfaJid^Xad^X^da, 

JM. 2 \Oai*\ 

P 6 =A p PV / ai -| 9 ai Xi [9 Q2 X 3 (ifr d ai ) (^ X 3 ) - (^1^)^X3^X3)]^ 

JR 2 Fai^-I 

+ V^ / i a Ql v^ ^^i^i^)^^^^) ~ a Q1 X 3 (i?id ai )(,9^X3)]<fa. 

JR 2 I C ai A I 

We control i-fc as before. Regarding P4 we use (17.71) to write it in the form P& = S\ + S 2 + S3 
where: 

5 1 = -A P PV [ T ^-^d ai X 1 (R 1 d ai ){d a2 X 1 d^X 1 )da, 

J R 2 \a ai x\ 6 

5 2 = -A p PV [ a ai 2 d ai X t (R t d ai ) (d a2 X 2 d±X 2 )da, 

Jm 2 \o ai x\ 

S 3 = -A P PV [ ^l^-d ai X 1 (R 1 d ai )(lo.t.)da. 

The identity (17. 5j) allows us to write P5 = S4 + 5s + S^ where: 

5 4 = A P PV / - Q1 2 d a2 X 1 (R 1 d ai )(d ai X l d^ 1 X 1 )da, 

5 5 = A p PV I .J Xl ^ 2 ]3 d a2 X 1 (R 1 d ai ){d ai X 2 dt 1 X 2 )da, 

JR 2 l^ai^l 

s 6 = V 3 ^ / _ ai ' a aa x 1 (JZia ai )(i.o.t.)cfa. 

Jr 2 l^a^i 

Next, we reorganize the sum in the form P4 + Pq = (S\ + S4) + (S 2 + S5) + (S 3 + Sq) where 
the term S 3 + Sq can be easily estimated. Regarding S\ + S4 we have 

Si + S 4 = A p PV [ .J* 1 * 2 d ai X 1 [d a2 X 1 (R 1 d ai )(di 1 X 1 ) - (Rrfa^d^X^X^da 

JR 2 |C ai A| 

r rfi X 

+ A P PV / la" 1 xrn dcvXlKRldaiWaxXid^X!) - d^X^RtfaJ^X^da 

Jm 2 IOc^AI 



and the commutator estimates gives us precise control. 
Let us consider now 



S 2 + S 5 = A p PV [ 1 J* 1 ^| , d Ql X 1 [d a2 X 2 (Ri d ai ) {8 A ai X 2 ) - (i? 1 a Q1 )(a Q2 X 2 ^ 1 X 2 )]da 

JR 2 |CaiA| 

r rfi X 
+ A P PV / 1 a" 1 y ?3 9a 2 X l [(R l d ai )(d ai X 2 di 1 X 2 ) - d ai X 2 (Rid ai )(dt 1 X 2 )]da 

Jm 2 \o ai x\ 

~ A p PV I T^^^ 1 A' 2 (i? 1 a ai )(^ 1 X 2 )da. 

7R 2 l^ajAI 
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Here again the commutator estimate control the first two integrals above, allowing us to 
conclude that 

2 = P 4 + P 5 + P e <P(\\\X\\\ 4 )-A p PV [ -J^—dt 1 X 2 (R 1 d ai )(di 1 X 2 )da. (7.14) 

Furthermore, inequalities (|7.13p . (|7,14p and (|7.12[) yield 

M l = O l + 2 + O z <P(\\\X\\\ A )-A p PV f —^-3 d 4 ai X ■ (RMidi^da, (7.15) 

and at this point we begin to recognize the Rayleigh- Taylor condition in the non-integrable 
terms. Let us return now to the term M 2 which can be written in the form 

[' B^ X N 

M 2 = A p PV / fli( ,% _.,, )di 1 {BR(X,oj) ■ d ai X))da, (7.16) 

JR2 \O ai A\° 

and whose most dangerous components are given by 

4vr J r2 I^X^ J r2 |A(a)- A(/3)| d 



_3A,.„„ A „ ,dtX-N 
with 



o a r fp- X N P 

5 = -^PV R 1 ( g y . 3 )(a) / S(a,/3)(X(a)-X(/3))A W (/3)-a ai X(a)da, 

(X(a)-X(/3))-(^ 1 X(a)-^ 2 X( / 3)) 



B(a,0) 



|X(a)-X03)|5 

''ai- 



O e = -^PV [ Rl &?LJL)( a )[ *ff *ff A ^(fl ■ d ai X(a)da, 

4vr J R 2 Ic^A^ 7 R 2 \X(a) - X{j3)\ 6 



p d x N 

O r = A p PV / R 1 {-^—- r )(a)d ai (BR(X,u J ) ■ d^X^da. 

JR 2 \0 ai A\ 

The remainder terms are less singular and can be estimated with the same methods used 
before. To deal with O4 we decompose it further O4 = P-j + P§: 



\ pv f R ( ^ X - N ) (n . ) f < X ^ ~ < X ^ 



p r = 7Z PV LM £ v| 3 )( g ) / Tv7 /_ //mi3 ^ M^(d ai X(f3)-d ai X(a))d(3da, 



P8 = ~rPy Rl{ IP y| 3 )(«) / |y, X y/ m ,3 " N {^d^^dPdXX, 

4vr J R 2 |<9 Ql X| d J R 2 |A(a)-A(« 



where in Ps we have used formula (|2.3|) to get uj A <9 Q1 A' = Nd ai Q. In the integral (with 
respect to (3) of Pj we have a kernel of degree —2 applied to 4 derivatives, which can be 
estimated easily. Next let us consider P§ = Qj + Q$ + Qg where 

Qr - ~^ PV L Ri{ ^xr ){a)d ^ x{a) • L — \xjaTx\w — d/WQ> 

\ /* /-)4 "V TV T {' 

Q 8 = -HPV R x { g y )(a) / ((a ai niV.^ 1 X)(a)-(a ai fiiV.^ 1 X)(/3))C(a^)d/3da, 

471" Jr2 |C Ql A|° J r2 
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and 

C(a ' ® = \X(a)-X(P)\ 3 ~ |VX(a)(a-/3)| 3 ' 
A f B^ X N 1 

In Qj we have 

,<9iX-iV\„ „„, „,. /■ N(a)d ai n(a)-N(P)d ai n(0) 






|X(a)-X(/3)|3 



tf/3 



giving us the appropriated control, which can be also obtained in Qg because the correspond- 
ing kernel has degree —2. Regarding Qg we have the expression 

* - t pv L R ^w )[ w^^' aN -^ x) - A( 8 "Zxf x >]da 

4vr J R 2 \d ai X\ 6 \o ai X\ 6 

Then we use (|9.2p to control the first integral above, and since A = R\d ai + R2d a2 (|9.4p we 
can also take care of the second term. 

With O5 one proceed as we did with Jq (|7.10p to get the desired estimate. 

Next we use (|2.3|) to catch the most singular terms in Oq which are given by 

ST = -£ PV I *(%^)<a) / ^M-^))AW). 8 .,X M 

8vr 2 7 R 2 |<9 ai X| d y R2 |X(a)-X(/3)| d 

a _ ^ pv f R ( < X - N ) ( ^ f (X(a)-X(f3))Ad Q2 XW).d ai X(a) 5 

510 = s^ pv L Rl (^W ){a) L |a»-aw ^n^U^- 

One may write 

4vr J R 2 \d ai X\ a J R2 \X(a) - X{(3)\ 6 

expressing the fact that we have a kernel of degree —1 applied to i9j 1 (9 Q , 2 J7 and, therefore, an 
integration by parts gives us the desired control as we did before. To treat 5§ we decompose 
further S% = T\ + Ti : 

a — r _ ,dtx-N 



where 



and 



T X = —±FV / R x { .g -3 )(a) / D(a,/3) • ^nG^X^da, 

4vr J R 2 |C Ql A| d J R2 

r (X(a)-X(/3)) VA»( a -/3) , 

T 2 = ^Py / ^i( ^ lX y| f )(a) N ^f ■ R 2 (d a2 nd 5 ai X)(a)da. 
An y R2 \o ai X\ 6 \d ai X(a)\ 6 
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In Ti we use the estimate for the operator (|9.8|) . The term T 2 reads as follows: 

A f B 4 X N N N eft X 

+ ^L PV / fll ( g . )[^_ . (i^X^fi^X) - (fl 2 3 ai )(,9 a2 • °* )]da 

The first integral above is easy to estimate, whether for the second one we use (|9.ip and (|9.4|) 
for the third. 

For the next term Sg one has Sg = T3 + T4 where 

/" (9^ A" TV 

T 4 = -^ / fl 1 (-^__)x>(^ 0)da, 

J R 2 \O ai Ji\° 

Proceeding as before we get bounds for T3 and the double layer potential estimates help us 
to control T4. 

For Siq one can adapt exactly the same approach used for Sg. Finally we have to deal 
with Oj which is given by 

f eft X N 

O r = -A^PV / BR(X,u) ■ d 4 ai X(R 1 d ai )(-^-— ir )da, 

Jr2 \o ai X\° 

after an integration by parts. Let us introduce the splitting O7 = Y2j k=i Uf where 

U* = -A»PV f BR j (X,u;)dt 1 X j (R l d ai ){ ^^ )da. 
Then the commutator estimates allows us to write Uj = VI- + lower order terms, where 
V} = -A li PV f BR j (X,u)d 4 ai X j -\ s (R 1 d ai )(d 4 ai X k )da. 

Using (|7.5p and (|7.7p one has 

iV 1 ^ 1 A 2 = iV 2 ^ 1 A 1 +l.o.t. 



so that V 2 becames 



V} = -A„PV 



BR 2 {X,uj)N 2 

\d ai x\* 



d 4 ai X 1 (R 1 d ai )(di 1 X 1 )da - A^PV [ f(R 1 d ai )(8 4 1 X 1 )da 

Jk 2 



where / is at the level of B^.X. Integration by parts in the last integral above allows us to 
conclude that 

V, 1 < -A,PV I BR " {X ^l N2 d i ai X 1 (R 1 d ai )(di 1 X 1 )da + P(\\\X\\U). 

JR 2 |C ai A|° 

With the help of (J73]) and ([721) we also S et 

iV 1 ^ 1 A 3 = iV3^ 1 A 1 +l.o.t. 
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and therefore 

Vi < -A,PV [ BR i {X ^l N3 dt 1 X l (R 1 d ai )(di 1 X 1 )da + P(\\\X\\U). 

JR 2 I Cai A | ° 

Using the two inequalities above we obtain 

Vi l + V} + Vi<-A»PV f BR ^^ N di 1 X l (R 1 d ai )(di 1 X 1 )da + P{\\\X\\\ A ). (7.17) 

Next let us observe that 

N 2 d A ai X x = N x d A ai X 2 + l.o.t., N 2 8 A ai X 3 = N z d A ai X 2 + l.o.t., 
which implies the estimate 

V? + Vi + Vi<-A»PV f BR ^\ Z N di 1 X 2 {R l d ai )(di 1 X 2 )da + P{\\\X\\\ A ). (7.18) 

Jm. 2 \0aiX\° 

Regarding V^ and V 2 the identities 

N^Xt = N x d A ai X z + l.o.t., Nad^Xs = N 2 d A ai X z + l.o.t., 
yield 

V? + VJ + Vi<-A»PV f BR ^] 3 N d A ai X 3 (R 1 d ai )(dt 1 X 3 )da + P(\\\X\\U). (7.19) 

Jr 2 \o ai x\° 

Finally (pUTD . (pU8D and (17TT9J) imply 

J2vf<-A„Pv[ BR{ *^] 3 N d A 1 X.(R 1 d ai )(d A 1 X)da + P(\\\X\\\ 4 ). 
jk=1 Jr 2 \o ai *\ 

Now we put together all those estimates ( (J7.16P - (J7.19P ) to conclude that 

M 2 < -A,PV f BR{ *^] 3 N d A ai X ■ (R 1 d ai )(d A ai X)da + P(\\\X\\\ 4 ), 
Jr 2 |C Q1 A|° 

and taking into account AT. 15f) we obtain 

L 3 = M X + M 2 < —Pvf ,, a v . 3 di 1 X-(R 1 d ai )(di 1 X)da + P(\\\X\\\ 4 ). (7.20) 

Finally we have to work with L5 which can be written in the following manner 

L 5 = L 5 - ]-PV [ 8 A ai X • 8 * 2 * A [R 2 (d A ai d a2 Qd ai X) - R 2 (d A ax d a2 ti)d ai X]da, 

2 Jr 2 \o a2 X\° 



where 



L 5 = \PV [ d A ai X- ^L^(R 2 d a2 )(d A ai Q)da 
2 Jr 2 \Oa 2 X\° 



Using the commutator estimate, once more, it remains only to consider L5, but let us point 
out that replacing the operator R\d ai by R 2 d a2 the term L3 (|7.1ip becomes L5. Therefore, 
proceeding exactly as we did before, one obtains inequality 

L 5 < ^—PV( jjr^d^X-iR.d^idi^da + PdWXWU). (7.21) 

M2+M1 Jr 2 lo^X^ 
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Introducing now the identity A = (Rid ai ) + (i?2^o 2 ) i n Q7.2UP an d (|7.21|) we get 

L 3 + L 5 < ]—FV l -^d^X.Kidi^da + PmXWU). 

H2+H1 Jr2 \d ai X\ 6 

Finally all the estimates so far obtained, beginning with (|7.9p . allow us to write 

\tM^\\1^)^-^T PV I T^^ d ^-Kd^X)da + P{\\\X\\U). (7.22) 



In a similar manner, using now equations ([2.9p . (|7.6p and (|7.8p instead of (|2.8p . (|7.5p and (|7.7p 
respectively, we obtain 

5|ll^lli»(*)< 7~ W -^di X-A(d* 2 X)d a + P(\\\X\\\ 4 ). (7.23) 

Being these two inequalities (|7.22p and (|7.23|) the main purpose of this section. 

8 Estimates for the evolution of ||F(X)||z,°° and R-T. 

In this section we analyze the evolution of the non-selfintersecting condition of the free surface 
as well as the Rayleigh- Taylor property, but in order to do that we shall need precise bounds 
for both V-Xt and Qf 

We shall estimate ||VX t ||^ by means of equality f|2.4[) to get 

\\VX t \\ Hk < P(\\X\\ 2 k+2 + \\F{X)f Laa + IHiVr 1 !!^), (8.1) 

for k > 2. In fact 

llVXiH^ < \\VBR(X,u)\\ H k + \\V{C l d ai X + C 2 d a2 X)\\ Hk 

and with the help of (16. ID we can handle both terms on the right. 
Next we shall consider the norms ||fit||#fc to obtain the inequality 

Untile < p{\\x\\1 +1 + \\F(x)f LOO + Hlivr 1 !!^), (8.2) 

for k > 3. To do that let us take a time derivative in the identity (j2.6|) to get 

n t {a,t) - A^V{n t )(a,t) = Aphfat) - 2A p d t X 3 (a,t), 

which yields 

Hfitllffi < C\\(I - A tl V)~ 1 \\ m (\\I 1 \\ H1 + \\d t X 3 \\ H i), 

and since we have control of \\(I — A ll T>)~ l \\ H i and [|9t.X3||#i it only remains to estimate 
||ii||#i. For that purpose let us consider the splitting I\ = J\ + J 2 + J3 where 

J 2 = ^ / (X(a)-X(a-/3))-(X t (a)-X t (a-/3)) , y /"? ~ y /" ~ f? . 5 -N(a-/3)n(a-l3)dp, 
47r J R 2 \X{a) - X(a - p)\° 

' h = ^ PV l \x(1~x( a ~al - N ^-^ a -^- 

2vr 7 R2 \X{a)-X(a- (3)1* 
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Proceeding as we did with the operator T2 (|9.6|) ( with X t instead of d a .Xk) one get 

\\Ji\\l* + Whhz < p(\\x\U + \\F(X)\\ L ~ + IWN^Wloo). 
Regarding J 3 we split further 

h = ^~ [ dp + ±- I dp = K 1+ K 2 . 
27r J\p\>i 27r J\p\<i 

m*)\ < nwii- / w-/)ijgf-wi ^, 

7|/3|>i 2vr|/^ 

Young's inequality yields 

II^H^ < \\F(X)\\ 2 L oo\\N t n\\v < C||F(X)||| 00 ||iV t || L2 ||Q|| L2 , 

and since we know that ||iVi||^2 < ||VX||£,cx)||VXt||^2, estimate ()8.1|) allows us to handle the 
terms K\. The estimate for K<i is similar to the one obtained for I2 (|9.13p in the appendix. 
Next we consider the most singular terms in d ai I\ which are given by 

t ! pv f d ai X t (a) - d ai X t (a - P) 

J A = —PV / — — — — — j3 N(a - P)Q{a - P)dp, 

2vr 7 R2 \X(u)- X(a- P)\ A 

J 5 = ^ / (X( a )-X(a-P))id ai X t (a)-d ai X t (a-P)) .^^~^ a ~^. 5 -N(a-PMa-P)dp, 
4vr J R 2 \X(a) - X(a- P)\° 

J 6 = 1 LpV [ 1 5°? ~ w" ~ £1 • 9 ai N t (a - PMa - P)dp. 
2vr J R 2 \X(a) - X(a- p)\ A 

because the remainder terms are easier to handle. Let us write J4 = K3 + K^ where 
ir ! dt/ / d ai X t (a) - d ai X t (a - f3) 

KA = ^ PV f ^(a)y(a-/3) . 

2vr J R 2 \X{a)-X(a-P)\ 6 

In Ks, the identity d ai Xt(a) — d ai Xt(a — j3) = f Vd ai Xt(a + (s — l)/3)ds ■ ft together with 
(|8.1|) gives us the desired control. Regarding K± we may observe its similarity with 73 (|9.7|) . 
so that an application to (|8.1j) yields the appropriated bound; J5 can be treated in a similar 
manner and Jg is analogous to J3. By symmetry, one could get the same estimate for d a2 I±, 
so that finally: 

\\n t \\ m < p(\\x\\l + \\F(x)f Loa + Ulivr 1 !!^). (8.3) 

Next, we will show how to deal with ||f^||#2. Using equation (|2.8p one gets 
d 2 ai tt t = -2A„d ai d t {BR(X,u) ■ d ai X) - 2A p d 2 ai d t X 3 , 



and with the help of (18. lj) . the last term above is properly controlled. To continue we shall 
consider the most singular remainder terms. Namely, in —d ai dt(BR(X,u}) ■ d ai X), we have: 

L 1 = -BR{X,u)-% 1 X t 

t l tdm f d ai X t (a) - d ai X t (a - P) 

L i = JZ PV — \Yu\ FT™ — «Fs — Aoj(a-/3)d/3-d Ql X(a) 

4vr J R 2 \X(a)-X(a- P)]* 
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L, = -^PV [ A(a,/?) * (a) y/" ^ Au;(a-{3)dP-d ai X(a) 

87r J R 2 \X{a) - X(a- p)\° 

where A(a,/3) = (X{a) - X(a - $)) ■ (d ai X t (a) - d ai X t (a - $)), 

L 4 = ±PV [ ff Q) ~ l\ a -ff A d ai u t {a - d)dp- • d ai X(a). 

Let us observe that ||Li||^2 < ||Bi?(X, w)||l°c .||9^ Xt||x,2, where both quantities have been 
appropriately controlled before. In L 2 and L3 we have kernels of degree —2, and therefore 
operators analogous to Tz (19. 7p acting on d ai X t . Therefore using (18. ip its control follows 
easily. In L4 we use the decomposition 

L 4 = -^-PV J dp + ^-PV J d/3 = M 1 + M 2 . 



27r J\p\>\ ' 2vr J\p\<x 
Thus an integration by parts yields 

||M 1 || i2 <C||F(X)|||o ||VX||| 00 ||w;t||i2. 



Formula (|2.3[) together with estimates (|8.ip and (|8.3p . provides the appropriated bound. 

Next let us expand (j2.3[) to obtain the most singular terms in M 2 which are given by the 
integrals: 

O l = -^PV I 3 a )~*( a ~»L A d a2 n{a - f3)d 2 ai X t (a - p)df> ■ d ai X(a), 
2vr J\fi\<i \ x («) - X \u-P)\ 

2 = -^PV f |5 (Q ?~5 (a "f?i3 A d ai d a2 n t (a - f3)d ai X(a - f3)d(3 • d ai X(a), 
2?r J\i3\<i\X{a) - X{a- j3)\ s 

3 = ^PV f ,5 (a ?~y/ a ~£| 3 A d a M» ~ P)d ai d a2 X t (a - (3)d(3 • d ai X(a), 
27r J\i3\<i \X[a) -X{a-P)\ 6 

Oa = ^PV [ ,y/ a ?~y/ a ~ff|3 A d 2 ai Q t (a - P)d a2 X{a - fldfi • d ai X(a). 

27r ./|/3|<i \X{a) -X{a-p)\ (i 

Estimate (|8.ip help us with the terms 0\ and O3, which can be treated with the same 
approach used for I 2 (I9.13P in the appendix. Let us write 2 as follows 

02 = ^1 1 w1"w a "fi A d ai daMa-P)(d ai X(a)-d Ql X(a-l3))dp ■ d ai X(a), 

27r J|/3|<i |^(a)-X(a-/3)|' :i 

which can be estimated integrating by parts in the variable /3\ using the following identity 

d ai d a2 n t (a-P) = -d^{d a2 Vt t {a-l3)). 

Let us point out that the kernel in the integral 2 has degree —1 and, therefore, one can use 
(|8.3p to control it. It remains to deal with O4 which is decomposed in the form O4 = P\ + P 2 , 
where 

Pi = i^Pvf |w (a ?~w a ~fi A d 2 ai n t (a-p-)(d a2 X(a-p-)-d a2 X( a ))dp- ■ d ai X(a), 

27r J\B\<i \X(a)-X{a-(J)\ i 



A MDT/ f X(a)-X(a-/3) a2 



27r J\0\<i\ x {<x)- x {<x-P)\ 
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Pi is estimated like O2. We rewrite P2 as follows 

-Vo„ f I X(a)-X(a-0) VX(a)-/3 



ft = -t PV ] w<1 \ \xia)-xi«-W ~ WxWW)^ a{a ~ m ' N{a) 



and this expression shows that the above integral can be estimated like T4 (|9,8 
Using (|8.3p we obtain 



IIC^II^ < P{\\x\\l + [|FWI|i- + lllivi" 1 !!^), 

and the identity 

d 2 a2 n t = -2A^d a2 d t (BR(x,uj) • d a2 x) - 2A p dl 2 d t x 3 , 

yields 

\\d 2 a2 n t \\ L2 < p(\\x\\l + \\f(x)\\ 2 loo + IHivr 1 !!^), 

that is: 

Untile < P(\\X\\l + \\F(X)\\l~ + HlJVr 1 ^-). (8.4) 

Next we consider third order derivatives 

d 3 ai n t = -2A^d 2 ai d t (BR(x,uj) • d ai x) - 2A p d 3 ai d t x 3 . 

Since (|8.ip gives us control of the last term, we will concentrate in the other one which is a 
much more diffecult character. In particular, for —d^dt(BR(X, u) ■ d ai X), the most singular 
component are given by 

L 5 = -BR(X,u)-% 1 X t 

t l pv f d 2 ai X t (a)-d 2 ai X t {a-^) 

Air J R 2 \X(a) - X(a-p)\ A 

L 7 = ^PV [ B(a, /?) * (a) ~ * (a Z® A U (a - P)d(3 • d ai X(a) 
8vr y R2 \X(a) -X(a-/3)| 5 

where S(a, /3) = (X(a) - X(a - (3)) ■ (d 2 ai X t (a) - d 2 ai X t (a - /3)), 

L * = k pv \ 1 ?ffi ~xt ~% A ^ (a - w • ^ x(a) - 

2vr J R2 I A (a) - A(a-/})| d 

Inequalities (|8.ip and (|8.4|) show how to handle Lj, % = 5, ...,8 as Lj, jf = 1, ...,4 respectively, 
then a similar approach for d^ fit allows us to get finally (|8.2p for k=3. The cases k > 3 are 
similar to deal with. 

Our next goal is to obtain estimates for the evolution of ||-F(a~)||loo and R-T. Regarding 
the quantity F(X) we have 

d |/3| (X(q, t) - X(a - /3, t)) • {X t (a, t) - X t (a - (3, t)) 

dt r^ )( a,p,t)- \X(a,t)-X(a-/3,t)\* (8.5) 

<(F(X)(a,(3,t)) 2 \\VX t \\ L ™(t). 

Then Sobolev inequalities in ||VA*t||£°o(£) together with (|8.ip yield 

±F(X)(a,M < F(X)(a,P,t)P(\\X\\l(t) + \\F(X)\\h.(t) + HM^Hl-W), 
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and an integration in time gives us 

F(X)(a,P,t + h) <F(X)(a,P,t)exp( j P(s)ds), 

for h > 0, where 

P(s) = P(\\X\\l(s) + \\F(X)f Loo (s) + |||Ar|-i|| LOO ( s )). 
Hence 

\\F(X)\\ L oo(t + h)< \\F(X)\\ L oo(t)ex P n P(s)ds). 

The inequality above applied to the limit: 

— F (A L - (*) = hm 

at h->o+ h 

allows us to get 

^\\F{X)\\ L ^{t) < \\F{X)\\ L ~(t)P(\\X\\l + \\F{X)f Loo + \\\N\- l \\ L ~). 

Next we search for an a priori estimate for the evolution of the infimum of the difference 
of the gradients of the pressure in the normal direction to the interface. Let us recall the 
formula 

a(a, t) = (p 2 - fj})BR(X, u))(a, t) ■ N(a, t) + (p 2 - p 1 )^^, t), 

to obtain 

d 1 _ cr t (a,t) 

dVo(a,ty~ a 2 (a,t) 
with at(a, t) = I\ + Ii where 

h = (Gu 2 - p})BR{X, u )(a, t) + (p 2 - p 1 )^, 0, 1)) • N t (a, t), 

h = (p 2 - p l )BR t {X, w)(a, t) ■ N(a, t). 



First we deal with ||ii||i,°° using the estimates (|8.1[) for VXt, and then we focus our attention 
on I2 using the splitting I2 = J\ + J2 + ^3 where 

1 „ Tr f X t (a) -X t (a-P) , 

J 1 = -—PV fv; V ^Awa-W, 

J2 = ^^y R2 ^(«)-^(«-/ 3 )) Aw («-/ ? ) |X(a)-X(a-/3)P ^ 

4vr 7 K 2 |X(a) -X(a-/3)| J 

The terms J\ and J2 are similar and can be treated with the same method. Let us consider 
j 1 =K 1 +K 2 + K 3 + K A where 

1 f X t (a)-X t (a-(3) 

Ki = --r- / 7^-^ t^t STT3 Aw(a-/3 )dp, 

^ J\p\>\ \X(a)-X{a- P)\* 
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1 f X t (a) - X t (a - B) , , , , 
K2 = j- ,*;; y , U -^ A w(a) - w(a - P ))d/3, 

4vr V|/3|<1 \ X ( a ) ~ X ( a ~ P)\ 

Ks = "J" / [ jy? S W M- |u Y /\ |3 ](X f (a)-X f (a-/3))Aa;(a)d/3, 

47r J|/3|<i |A(a)-A(a-y3)| d |VA(a)-/3| d 

1 ™, f X t (a) - X t (a - P) , > ,„ 



First we have 

||tfi||i» < C||F(X)||| 00 ||VX t || L o || w || L2 ( / 1/31-^/3) 1/2 

■^|/3|>l 
giving us an appropriated control. Next, we get 



\\K 2 \\ L °°<C\\F{X)\\l~\\VX t \\ L oo\\Vu\\ L o° I \f3\~ l dp, 

J\P\<i 

and an analogous estimate for K3. Therefore, Sobolev's embedding help us to obtain the 
desired control. Regarding K^ we have 

**J m< i |VX(a)-/3| 3 l ' P 

Inequality (19.151) yields 

\\K4l°° < C||VX||| 00 |||iVr 1 ||| 00 || w || LO o||VX t || c , I \p\- 2+5 d[3, 

■/|/8|<1 

and the control || VXt|| C 6 follows again by (|8.1|) and Sobolev's embedding. Next let us continue 
with J3 = A'5 + Kq where 

k 5 = =±pv[ 2^t£j?L^) A(d Pl ((nd a2 X) t (a-P))-d P2 ((nd ai X) t (a-P)))dp, 

47r J\p\>i \X(a)-X{a-P)\ 6 

K « = -T PV I i^",?"^ ^-^ 

47r V|/3|<i \X[a) -X{a-p)\ A 
Integration by parts yields 

II^IIl 00 < c||^ ? (^)||z / oo||vx||ioo(||j7||io ||vx i ||ioo + 1 1 r2t 1 1 x,°° 1 1 v j^t 1 1 z,°° ) , 

where 4-kC = fiai >1 \P\~ S dp + Jig\ =1 dl(fi), and we may use (|8.2p to estimate ||fii||z,°°. With 
Kq we introduce a similar splitting to obtain 

\\K 6 \\ L «> < P(\\X - (a,0)\\ C 2 + \\F(X)\\ L ~ + WiVriL-OIMIc*. 

Then it remains to estimate ||wt||c7*5 f° r w bich purpose we use formula (J2.3H and inequalities 
(|8.ip()8.2p . Therefore we have the estimate: 

|(^-r) < ^-rPdl^lUW + [|FW||l-(*) + \\\N\- l \\ L ~(t)), 
at a{a,t) a z {a,t) 

and proceeding similarly as we did for F(X) we get finally: 

^Ik-'IU-W < llOli-CWII^IUW + ||FWIlL-(t) + IH^r'lUccCt)). 
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9 Appendix 

Here we prove first some helpful inequalities regarding commutators of Riesz transform (Rj, 
j = 1,2) with several differential operators. Next we analyze the singular integral operators 
associated to the non-selfintersecting surface which appears throughout the paper. But the 
main goal of this section, however, is to simplify the presentation of the main result. 

Lemma 9.1 Consider/ G L 2 (R 2 ), and g G C^ S (R 2 ) withO < S < 1. Then for any k, I = 1,2 
we have the following estimate 

\\(R k d ai )(gf)-g(R k d ai )(f)\\ L , < C\\g\\ c n.s\\fh 2 - (9.1) 

An application of the above inequalities to the operator A = (R\d ai ) + (i?2<9 Q2 ) yields 

\\A(gf)-gA(f)\\ L 2<C\\g\\ cl ,s\\fh*- (9-2) 

For vector fields we have 

Lemma 9.2 Consider /, g : R 2 -»■ R 3 vector fields where f G L 2 (R 2 ) and g G C 1 ' 5 (R 2 ) with 
< 5 < 1. Then for any k,l = 1, 2 the following inequality holds 

(gAf)-(R k d at )(f)da\<C\\g\\ cl , s \\f\\h- (9-3) 

Proof: Denoting with / the integral above and since the operator R k d ai is self-adjoint we 
may write 



I = / h[(Rkd ai ){92h) - 92(Rkd ai ){h)}da + / f 2 [(Rkd ai )(g 3 fi) - g 3 (R k d ai )(h)]da 

+ / h[{Rkd ai ){gih) - gx{Rkd ai ){h)]doc. 

Then estimate flEE]) yields ([Ojl . 



Lemma 9.3 Consider f G L 2 (R 2 ) and g G C 1 ' 5 (IR 2 ) with < S < 1. T/ien /or any j, fc, / 

1,2 the following inequality holds 



R 3 {f){R k d ai ){gf)da\ < C|| 5 || cl ,,||/||| 2 . (9.4) 

Proof: Let J be the integral to be bounded, then we have 
J= [ RjiMRkdaJigf) - g(R k d ai )(f)]da - I [R j (fg)-gR j (f)](R k d ai )(J)da 

+ / R 3 {fg){R k d ai )(f)da 

Jm. 2 



Since R* = —Rj and R k d ai is self-adjoint we get 



J= \ I Rj(f)[(Rkd ai )(gf)-g(R k d ai )(f)]da-l- [ [R j (fg)-gR j (f)](R k d ai )(f)da. 

2 Jr2 2 J R 2 

An integration by parts in the second integral above yields 

J = \ f Rj(f)[(Rkd ai )(gf) - g{R k d ai ){f)\da + \ f [{R d ai ){fg) - g(R 3 d ai )(f)](R k )(f)da 



{d ai g)R 3 (f)R k {f)da, 



allowing us to conclude the proof. 
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Lemma 9.4 Let us define for any j = 1, 2 and k = 1, 2, 3 the following operators: 

T x {d aj f){a) = PV I ^(«)-^(«-^) g /(a _ p )d p (9 . 5) 

Jr2 |A(a) - X(a- fj)^ 

f d a .X k (a)-d a .X k (a- p) , 
T 2 (f)(a) = PV a > k a k{ P, f(a-(JW, (9.6) 

Jm? \X(a) - X(a- P)\ 6 

r 3 (/)(a) = J>W i/ff'y^'oL ^ ( 9 - 7 ) 

Jr2 l^(a) - A(a-/J)| d 

-r^ fV ^ pv /" f (*(<») -*(£)) VX(a)-(a-/3) x a 

T 4 (5 Qj /)(a) = PV^ ( |x(a) _ x(/3)|3 " | V X(a)-(a-/3)|3 )^ /(/3) ^ a ' (9 " 8) 

iu/iere VI(a) • /3 = d ai X{a)Pi + d a2 X(a)P 2 . Assume that X(a) - (a,0) E C 2 ' 5 (M 2 ), and 
that both F(X) and \N\~ l are in L°° where 

F{X){a,P) = \P\/\X(a)-X(a-/3)\ and N(a) = d ai X(a) A d a2 X{a). 

Then the following estimates hold: 

\\Ti{d a .f)\\ L * < P(\\X - (a,0)|| cM + \\F(X)\\ L ~ + |||iV|- 1 ||L-)(||/||L 2 + \\d a J\\ L ^ (9-9) 

||r 2 (/)|| L2 < P{\\x - (a,o)|| c2> , + \\f{x)\\ loo + IHivr 1 !!^)!!/!!^, (9.10) 

nr 3 (/)ii L2 < p(\\x - (a,o)\\ c ^ + \\nx)u«> + mr'h^wfWm, (9.11) 

||T 4 (/)|| L2 < p(\\x - (a,o)|| c2 ,, + \\F(x)\\ Leo + lllivr 1 !!^)!!/!!^, (9.12) 

with P a polynomial function. 

Proof: To estimate the first set of operators we consider first the splitting 

Ti(d a .f) = PV i dp + PV [ dp = h + I 2 

J\p\>i J\P\<i 

and an integration by parts allows us to write I\ = J\ + J 2 + J3 where 

f -d a .X k (a- P) 

J\/3\>i \X(a) -X(a- P)\ 6 

T . /• (X fc (a) - X k (a - P))(X(a) - X(a - p)) • d aj X(a - P) 

J 2 = 3 / r^-r^ —, ^ /(a - P)dp, 

J\fi\>i \X{a)-X{a- P)\ b 

and 

/• X fc (a)-X fc (q-/3) ^ 

J\B\=i \X{a) -X(a- p)\ A 



(9.13) 



The above decomposition shows that 

|/ 1 |<C||X-(a,0)|| cl ||F(X)||| oo ( " K "~ 3 ™ dp+ I/O 



|3 ,/■ |/(«-/3)| ^ 

'|/3|>1 I/ 5 ! 3 ' 7|/3| = 1 



and then Minkowski's inequality gives the desired control. 
Regarding L 2 we write I 2 = J4 + J5 + J@ with 

, _ f X k (a)-X k (a-P)-VX k (a).p 

'U- I 1 w x 777 3TT3 O aj J{a- p)dp, 

J\fi\<i \X{a)-X{a- P)\ 6 
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J 5 = VX k (a) ■ [ 0[ I —3 - ——L—-]d a .f(a - P)d/3, 

J\p\<i \ x { a ) - X (a-P)r |VA(a)-/3| d 

J 6 = VX k {a) ■ PV [ 

J\0\<1 

It is easy to see that 



\xjy( \ R\3 d <*jf( a ~ P) d P- 



J*<\\X-(a,0)\\ c i,s\\F(X)f Loa / "' L* d/3, (9.14) 



I \d aj f(a -P)\ 

|/3|<1 W 

and therefore that term can be estimated also with the use of Minkowski's inequality. 
Some elementary algebraic manipulations allows us to get 

J 5 <C\\X-(a,0)\\ln, [ [(F(X)(a,/3)) 4 + l/3 ' 4 ] lgSj^M djg| 

J\fi\<i \vA{a) ■ p| \p\ 

and then the inequality 

^br 2||vx|WI|Arl|ll ~ (915) 



yields for J5 the same estimate (19.14|) . 
The term Jg can be written as 



J 6 = VX k {a) ■ PV [ ^fl-d aj f( a - [3)d/3, 

J\P\<1 



|/3|2 ° a 



where 

(t) S(a, A/3) = E(a, /3), VA > 0, (w) E(a, -/3) = -E(a, /3), 

and 

(in) sup|E(a,/3)| < 8[|VJT][|,oo |||iV| -1 |||,oo, 

as a consequence of (|9.15p . 

Here we have a singular integral operator with odd kernel (see [8] and [T7]) and therefore 
a bounded linear map on L 2 (M 2 ) giving us 

l|j 6 || L2 < c||vx||l^|||Ar1!oc||<9 aj /|| L2 . 

For the family of operators 7^(7) (a) we use the splitting 72(/) = X3 + Ii where 

d aj X k (a) - d aj X k (a - /3) 
w>1 \X(a)-X(a-/3W 



^ = I a W,J v l. m,3 - /(" - PW- 



Easily we get 



/ 3 <2||X-(«,0)|| cl ||F(X)||| oo / l/( " ' R / )l d0, 

J\0\>1 \P\ 



while for I 4 we proceed with the same method used with I 2 , replacing now X k (a) by d a . X k (a) 
and d aj f(a - (3) by f(a - /3). 

Next we shall show that the operator T3 behaves like A = (—A) 2. To do that we split it 

as 1$ + Iq where 

f(a)-f(a-(3) 



0>1 \X(a)-X(a-P)\ 3 
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dP, 



can be easily estimated by 



h < \\F(X)f Loa (2^f(a)\ + [ l/( " 3 /3)l d/3). 

J\0\>i \P\ 



I/3|>1 
The other term is written in the form Iq = Jj + Js where 

■V|<i l^(a) - X[a - P)\ A |VX(q) • fi\ A 
The identity 

/(a) - /(a - /3) = /3 • / V/(a + (s - l)0)ds 

JO 

allows us to treat J7 as we did with J5. To estimate Jg the equality 

1 Q / h \ a ( fa 



will be very useful. After a careful integration by parts it yields 

J\8\<1 |VX(a)-/3| d J m=l |VX(a)-/3| d 

The principal value in J$ is treated with the same method used for Jq and since the integral 
on the circle is inoffensive, so long as |iV| -1 is in L°°, the estimate for T3 follows. 
For the remaining operator one integrates by parts to get T4 = I7 + Is where 



I 7 = PV Pi(a,(3)f(a - (3)dp, h = PV P 2 {a^)f{a - (3)d(3 
Jk 2 Jk 2 

with 

d aj X(a) d aj X(a-P) 

Pi{a,/3) 



\VX(a)-f3\ 3 |X(a)-X(a-/3)| 3 
and 

\X{a) - X(a - 0))(X(a) - X(a - &)) ■ d a .X(a - 0) 



P 2 (q,/3)=3- 



\X(a)-X(a-P)\ 5 
yX(a)-P((VX(a)-P)-d Qj X(a)) 

|VX(q)-/3| 5 



Next we will show how to treat I7, because the estimate for ig follows similarly. In Pi we 
introduce the further decomposition: P\ = Qi + Q2 where 

n R Y( u 1 1 1 n _ d aj X{a)-d aj X(a-P) 

Wi a "i*W[\ VX ( a y p\s \X{a)-X{a-p)\^ ^ 2 \X(a) - X(a - /3)| 3 " 

And since the kernel Q2 has already appeared in the operator 71, it only remains to control 
Jg which is given by 

J 9 = d a .X(a)PV [ Qi (a, /3)/(a - /3)d/3. 



The following decomposition 

J 9 = d a .X(a) I dp + d a .X(a)PV I df3 = K x + K 2 

J\0\>i J\0\<\ 



'I/3|>1 J\P\<1 
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shows that the term K\ trivializes. Regarding K2 let us write 

{\A\ A + \B\ 2 \A\ 2 + \B\ A ){A + B)-{A-B) 

|A| 3 |B| 3 (|A| 3 + |B| 3 ) 

where 

A(a, 0) = X(a) - X(a - (3), B(a, p) = VX(a) ■ /3. 

This formula shows that inside Q\ lies a kernel of degree —2. Then let us take Qi = S\ + £2 

where 

_ 3\B\ 4 B ■ (A - B) _ 3B ■ (A - B) 

2 ~ \B\9 ~ \Bj 5 ' 

Next we check that the kernel Si has degree —1, and therefore is easy to handle. Finally we 
have to consider the kernel S2 appearing in the integral L 



(VX(a) ■ p) ■ (X(a) - X(a - (3) - VX(a) ■ (3) 
W\<i |VX(a)-/3|' 



L = 3d ai X(a)PV J w^^; hj ^j-^-^-v^, hi f{a _ p) ^ 



To do that we introduce a further decomposition L = Mi + M2, with 

vi /■ (VX(a)-/3)-(X(a)-X(a-/3)-VX(q)-/3-^-V 2 .. 
Mi = 3d aj X(a) I — , nv/ ^ fl | S -f{a-(3)d(3 

and 



2 7|/3|<i \yX{a)-p\° 

where i/3 • V 2 X(a) • /3 is the second order term in the Taylor expansion of X. It is now easy 
to check that 



MiKCWVXWl^WX - {aM\c^\\W\- X \\U / U LoY l d(3 



|/9|<i I0I 2 "' 

Then we also check that M2 is controlled like Jq throughout the estimate 
\\Mo\\t2 < CllVXllfoollV 2 ^! 1 - in ivrl— 1||4 



which allows us to finish the proof. 

Remark 9.5 Having obtained the a priori bounds of the precedent sections, we are in position 
to implement successfully the same approximation scheme developed in f^j to conclude local 
existence. 
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